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Abstract 



We investigate qualitative properties of local solutions u{t. x) > to the fast diffusion equation, 
dtu = A(ii/^)/m with m < 1, corresponding to general nonnegative initial data. Our main results 
^ are quantitative positivity and boundedness estimates for locally defined solutions in domains of 

the form [0, T] x WK They combine into forms of new Harnack inequalities that are typical of 
fast diffusion equations. Such results are new for low m in the so-called very fast diffusion range, 
. precisely for all m < rric = (d — 2) /d. The boundedness statements are true even for m < 0, while 

the positivity ones cannot be true in that range. 
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Introduction 



We study qualitative and quantitative properties of solutions u = u{t, x) of the nonlinear diffusion 
equation 



in the whole parameter range —oo < m < 1, where it is called Fast Diffusion Equation (FDE). We 
consider local nonnegative weak solutions, defined in an open cylinder Q of space-time M x M'^ with 
d > 1. Note that the factor 1/m'm the last expression is inessential when m > (up to a time rescaling, 
t' = t/ni) but becomes essential for m < 0, in order to obtain a parabolic equation; for m = the last 
expression has to be written as dtu = Alog(ti)[^ 

Assuming the basic existence and uniqueness theory, [Ij], [31], we are interested in the qualitative 
properties of the solutions such as boundedness, positivity, and Harnack inequalities. For the FDE these 
properties depart from the properties of the linear Heat Equation (case m = 1), [3l], and even more 
from the Porous Medium Equation (case m > 1) , [35]. Moreover, they are still partially understood 
when m is far from 1, precisely for m < rric where rric = {d— 2)/d is called the first critical fast diffusion 
exponent. Our goal here is to obtain bounds from above and below for the solutions in that low range 
of exponents. We look for precise quantitative versions based on local estimates. Such estimates should 
be of interest in developing a general theory of this equation in the detail that is already known both 
for m > 1 and for nic < m < 1. 

Precedents and problems 

The existence and uniqueness of weak solutions of the initial value problem and other standard initial 
and boundary value problems for the FDE, as well as the main qualitative properties of the solutions 
(such as the ones already mentioned, or the asymptotic behaviour), are by now well understood when 
m is close to one, more precisely in the so-called good parameter range: rric < m < ij^ To be specific, 
when the problem is posed in the whole space, weak solutions are uniquely determined by their initial 
data if uq is a locally integrable nonnegative function, or even a locally finite Radon measure. In that 
case, the solution is C°° smooth and positive for all x £ and t > 0, and the initial data are taken 
in the sense of initial trace, [23|, [27^, [Hj. Solutions are bounded for data uq £ LP{W^) for any p > 1, 
and even for data in the Marcinkiewicz spaces MP(M'^), p > 1, [31]. They are locally bounded under 
the very mild restriction that uq is Radon measure, even if it is not globally finite. 

The theory of the FDE has been much less studied until recently in the subcritical fast-diffusion range 
m < rric, even under the condition m > 0, since essential difficulties have been found in the different 
chapters of the theory, like existence, uniqueness, and regularity. Note that < m < rric is possible 
only if d > 2. We refer for background to the book [31] that discusses in some detail the range m < rric, 
even for m < 0, along with the cases m > rric. Let us give an idea of the difficulties that arise and that 
we address in our work below: 

Boundedness. Though weak solutions with data in the spaces LP(M'^), 1 < p < oo, exist and are 
unique for < m < 1, counterexamples show that for m < rric these weak solutions need not be 

^We will always interpret u"^ /m as log(«) when m — Q. In the whole paper, V indicates the gradient operator, V- the 
divergence operator, and A the Laplacian operator, all of them taken with respect to the space variables, x £ . 
■^With the extra restriction m > if d = 1, the case — 1 < m < and d = 1 being somewhat different, cf. [31] . 
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bounded, and as a consequence they are not smooth. The simplest such example seems to be the 
separate-variables function 

^(^'-;^'^o) = c^-^^^^ (0.2) 

For every m < nic, even m < 0, there exists a suitable constant c(m, d) > such that [/ is a weak 
solution of the FDE in the cylinder Q = {0,T) x M'^, cf. [31], page 80, but obviously the solution 
never improves its initial regularity until it extinguishes in finite time. The precise space regularity is 
U{-,t) G hf^^OR.'^) for all p < pc, where the critical integrability exponent is pc = d{l — m)/2, which is 
larger than 1 precisely for m < rric, i.e., in the subcritical range. 

There is a positive result concerning boundedness, that is also tied to the exponent pc '■ solutions with 
initial data in LP(M'^) with p > pc become bounded and C°° smooth for all positive times as long as 
the solution does not disappear. This smoothing effect happens for all p > 1 if m > rric, for p > 1 if 
m = rric (in the last cases there is no problem of disappearance). The results are sharp, cf. [31^. 

Extinction in finite time, EFT. The above example exhibits another typical feature of the Cauchy 
problem for m < rric, namely, the possible lack of positivity due to EFT. The occurrence of EFT depends 
on the type of problem we consider. 

In the case of the Cauchy problem posed in with d > 3, Benilan and Crandall gave in [3] a proof of 
the extinction in finite time, EFT, of solutions of the FDE in the range < m < rric when uq G L*'(]R'^) 
with p = Pc- It is proved in [31] that EFT occurs for the solutions with m < rric for all functions with 
initial data in the Marcinkiewicz space MP=(M'^), hence in LP=(M'^). We recall the EFT does not happen 
for the Cauchy Problem when m > rric- 

In the case of the Cauchy-Dirichlet problem posed in a bounded domain with zero boundary data, 
EFT happens for all < m < 1. There is an interesting functional connection: we can show that EFT 
occurs if we have a global Poincare and a Sobolev inequality, and this result can be extended to more 
general settings, such as Riemannian manifolds, as it has been done by the authors in [7 . On the other 
hand, Benilan and Crandall's proof for the Cauchy problem is based only on the Sobolev inequality, 
but it holds only in the lower range m < rric- 

Harnack inequalities. Concerning finer regularity properties, the possible occurrence of EFT 
is compatible with the fact that nonnegative bounded solutions are positive, and consequently C°° 
smooth, as long as they are not identically zero, i.e., before extinction. However, the existence of 
EFT for low m is tied to the breakdown of the standard forms of Harnack inequalities, which are a 
strong tool in developing a regularity theory. Obtaining some kind of Harnack inequality is therefore 
a main research issue for m < rric and has been an open problem for some years. More specifically, we 
concentrate on parabolic lower Harnack inequalities of the type called Aronson-Caffarelli estimates [1] , 
and examine their consequences to obtain quantitative forms of positivity. An extension work has been 
done in ^ for rric < m < 1 but the method collapses for m < rric due to the very different properties 
of the solutions. As a consequence of our local smoothing effect and of positivity estimates, we will 
obtain some intrinsic Harnack inequalities of forward, elliptic or backward type, which are new in this 
range. 

In a recent preprint |18j, DiBenedetto, Gianazza and Vespri study the validity of intrinsic Harnack 
inequalities in the good range m > rric and show, with an explicit counterexample, that any kind of 
Harnack Inequality, intrinsic, elliptic, backward and forward can not hold if m < rric, for a fixed size 
of the intrinsic cylinder, that is, if we fix the size of the parabolic cylinder "a priori" in terms of the 
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value of u at the center of the cyHnder {to, xq) . They leave as an open problem to find which kind of 
Harnack Inequalities, if any, are typical of the very fast diffusion range < m < me- In this paper we 
give an answer to this intriguing problem. 

Very singular range. Most the literature has avoided the cases m < 0, where the diffusivity 
D{u) = is very singular at n = 0. Recently, it has been shown that a large part of the theory of 

the subcritical range goes over to this very singular range, on the condition of working with solutions 
that "are not too small". See |30j among the older references, then [13], and the books [13], j31j for a 
more complete reference. Note that this recovers a subcritical range for dimensions d= 1,2, and also 
that we can study the interesting log-diffusion problems where m = 0, cf. [23], [33: and the references. 

More specifically, there is an extension of the results called smoothing effects, whereby data in L''(]R'^) 
with p > pc imply bounded solutions for all t > 0, and also the extinction in finite time for data in 
M^(M'^), p = Pc- But a very different situation happens for data in LP(]R'^) with 1 < p < Pc, which is 
called immediate extinction, whereby the solutions obtained as limit of any reasonable approximation 
are identically zero for all t > 0. This makes it difficult to think of a general study of positivity. 
Immediate extinction happens for the Cauchy-Dirichlet problem posed in a bounded domain with zero 
boundary data for all m < 0, d > 1. Our study of this range is confined therefore to upper estimates. 

Comparison with elliptig problems. Part of the difficulties of the FDE in the lower range of m can 
be explained by the intimate relation of the equation with semilinear elliptic theory. This remarkable 



connection will be briefly explained in Subsection 4.3 



Results and organization 

Our work focuses on a peculiar feature of the FDE, which is the existence of very strong local estimates. 
This was presumably first mentioned in the paper by Herrero and Pierre [23], 1985, who get solutions 
in the whole range < m < 1 under the sole condition on the initial data uq G Lj^^{W^). Much of 
the subsequent work has been influenced by the local character of the equation. Here, we want push 
this idea to its final consequence concerning two different areas: the question boundedness of local 
solutions, and the question of positivity of nonnegative solutions, measured quantitatively by so-called 
lower Harnack inequalities. We will then combine the local upper and lower estimates, into a full form 
of Harnack inequality. While the boundedness results hold for all m < 1, positivity estimates are 
confined to < m < 1 because of the possible occurrence of immediate extinction. As we have said, 
the main interest of our results lies in their application in the subcritical range, m < m-c. They are also 
new for the critical exponent m = nic- 

Let us be more specific about the contents of the paper. It is divided into three main parts. 

(I) The study of positivity and lower Harnack inequalities, both of local and global type. The first 
main contribution of the paper is a parabolic lower Harnack inequality of the Aronson-Caffarelli type 
that is presented in Section[T| along with a detailed comparison with the forms available for other ranges 
of m. We devote Subsection |1 . 1 1 to prove the lower estimate, Theorem |l.l[ for a minimal problem. This 
is extended in Subsection |1 .2| to general solutions. We then show that in the range m-c < m < 1 we can 
further eliminate the presence of the extinction time and recover stronger estimates that are known in 



that range. Subsection 1.4 discusses upper bounds for the extinction time T in terms of norms of 
the data, which give an alternative type of lower bound in the range where estimates depending only 
on L^-norms of the data are not true. 
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(II) The study of local upper bounds. This takes two forms: the first is the control of the evolution in 
time of some spatial Lj^^ norms, which is performed in Section 2.1 Then, we get a local in space-time 
version of the smoothing effect from Lj^^ into L^^, an important regularity result that opens the road 
to higher regularity and was known for m > rUc, and is false in general for m < rric. We show in this 
paper that the estimate holds m < rric, on the condition that p must be large enough. We finally obtain 
the finest local upper estimates, called local smoothing effects, in the form given in Theorem |2.1| , just 
by combining the space-time smoothing effect and the L|^^ obtained in the first Section 



2.1 



(III) Parabolic Harnack Inequalities. In Section [sj we combine the local upper and lower estimates 
obtained in Parts I and II in the form of parabolic Harnack inequalities of forward, backward and 
elliptic type, together with an alternative form. 

To conclude, we sketch a panorama of the obtained local estimates depending on the ranges of m, 
together with general remarks, some related open problem and a short review on related works. A final 
Appendix contains some useful technical results. 

Notations. We will work with weak solutions u > of the FDE with m < 1, defined in a cylinder 
Q = J7 X (To,Ti) for some domain O C M'^ and Tq < Ti. Usually, we take To = 0,Ti = T. Ti can be 
infinite and 0, can be the whole space. In view of existing theory we may assume that the solutions 
are positive and smooth as long as they do not extinguish identically. We will be mostly interested in 
the local theory where the space domain is bounded and the boundary conditions are not taken into 
account. In deriving local estimates it will be often sufficient to take as space domain a ball, which we 
will denote hy B = Br{xq) or B = B\fi(xo) for some A > 1. We will frequently consider the annulus 
region A^^x = Bxji \ Br. As indicated before, we put 

d-2 d(l-m) 
^ d ' ^ 2 

We have pointed out that pc > 1 if and only if m < rric. We will take integrability exponents p > 1 if 
m > rric, p > Pc ^ w-c- Moreover, for p ^ pc we set 



= i ^ , (0.3) 

^ 2p - (i(l - m) ' ^ ^ 



which is positive if and only if p > pc. 



1 Part I. Local lower bounds 

The first part of the paper addresses the question of quantitative estimates of positivity. The exponent 
range in this part is < m < 1, since it is well known that the FDE does not admit solutions of the 
Dirichlet problem with zero boundary data when m < 0, thus blocking any possibility of a general local 
positivity theory in that range \30\ |3T] . Our main contribution is a parabolic inequality in the spirit of 
the one obtained by Aronson and Caffarelli [T] in their path-breaking paper for m > 1, and the ones 
produced by the authors in [8] for rUc < m < 1. The purpose of such formulas is giving quantitative 
information on the positivity of solutions at later times in terms of information on L-P norms of n at a 
former time that we take as t = 0. This is why they are called parabolic lower Harnack formulas. 

We take < m < 1 and consider a n be a local, nonnegative weak solution of the FDE defined in a 
cylinder Q = (0,T) x J7, taking initial data u(0,x) = uq{x) in 0, and having finite extinction time T. 
We make no assumption on the boundary condition (apart from nonnegativity) . For ease of proof we 
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will assume that the solutions are smooth so that the different computations and comparison results 
are valid. This assumption is then eliminated by approximation, which is justified according to known 
theory. 



Theorem 1.1 Let < m < 1 and let u he the solution to the FDE under the above assumptions. Let 
xq he a point in i7 and let d{xQ, dO,) > 3-R. Then the following inequality holds for all < t < T 

R-'^ I uo{x)dx<CiR-'^/^^-"'h^+C2Tr^R-^t-^u"'{t,xo). (1.1) 

JBr{xo) 

with Ci and C2 given positive constants depending only on d. This implies that there exists a time t^, 
such that for all t G (0, t*] 

u'^{t,xo)>C[R^~^\\u4,^.^Bn)T~^ t^- (1-2) 
where C[ > depends only on d; depends on R and \\uq{x)\\i^i(^Bb.) ™^ on T. 



Simplified version. The dependence on the parameters makes the formula apparently complicated. 
But it can be reduced to a simpler, equivalent one. Actually, we may assume that xq = by translation. 
Given i? > and M = '^oi^) dx > 0, we use the rescaling 

n{t, x) = ^ u 1) , r = R2-<i(i-m)M^-m^ (13^ 

to pass from a solution with mass M in the ball of radius i? to a solution u with mass 1 in the ball of 
radius 1. So we only need to prove the version with M = i? = 1 to get the full version. The scaling is 
simpler for m = where r = M^~^. Of course, the extinction time has to be rescaled accordingly, 

T = ij2-d(l-m)^l-mji_ 



Improvements. As stated, estimate (1.1) applies only to solutions with finite extinction time, and 
it involves the value of the extinction time T in an explicit way; both things can make it impractical. 
However, a simple comparison argument shows that we only need to estimate from below any subso- 
lution. In particular, we may replace the solution under consideration by the solution of the problem 
with initial data uq{x)xbii(xo){^) ■• ^^"^ 'zeio Dirichlet boundary conditions on x G OB^ji^xq). Let us call 
this problem minimal problem for the given data. The extinction time of the corresponding solution 
will be called the minimal life time of such domain and data, Tm{uo, B). Clearly, Tm{uQ, B) < T{u). 



Corollary 1.2 The above positivity result holds with T (u) replaced by the minimal life time Tm^UQ, B), 
u is defined in Qt, and the estimate applies for < t < T' with T' = min{T, Tm}- 



This modified result is specially interesting in the range 1 > m > rric where the solutions of the 
Cauchy Problem do not vanish. On the other hand, it is known that is finite if uq satisfies some 
local integrability conditions I32|. 
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Comparison with the estimate for the PME and other FDE 



The PME. Let us write Aronson-CafFarelli's result p| for m > 1 with a similar notation: 

JBr(xo) 

We recall that this formula is valid for all nonnegative weak solutions of the PME defined in the whole 



space. The form of the first term in the right-hand side is the same in both results, (1.1) and (1.4). 
This term plays the role of blocking the positivity information when it is large relative to the left-hand 
side integral, and allowing for such information when it is small. The critical time at which we begin 
to get positivity information is obtained by making this term a fraction of the left-hand side, i.e., for 

= c(m,d)||no||[7('S,(.o))^'+'^'"-'^- (1-5) 

But since the exponents have just the opposite sign in the above expressions for m > 1 and m < 1, 
the consequences are qualitatively very different: the information on positivity happens for us when 
t is smaller than t^:, while for the PME it happens when t is larger. This is in accord with the basic 
properties of these equations, which the present inequalities faithfully reproduce. Rescaling allows to 
check the inequality only at t = 1 for R = 1, and in that case we only have to prove that there are 
constants Mq = Mo(n, m) and k = k{m, d) such that for M ^ Mq 

^(0,1) ^ A;M2/('^(™-i)+2). (1.6) 

As to the second term, it is different. We cannot expect to have the A-C term in the range m < rric 
since then the exponent of u would be negative. In fact, the proof of [1] uses conservation of mass that 
is not valid for the fast diffusion equation in the low m range. 

The good FDE. The validity of the Aronson-Caffarelli formula was extended by the authors in [8j 
to local solutions of the FDE in the good exponent range rric < m < 1, and the already mentioned 
sign change in the exponents implies that we get good lower estimates for < t < t*. Moreover, we 
can continue these estimates thanks to the fortunate circumstance that we have further differential 
inequalities, like dtu > —Cu/t in the case of the Cauchy problem, which allow for a continuation of 
the lower bounds for t > t^, with optimal decay rates in time. The final form is 

u{t, x) > Mr{xo) H{t/tc), Mr{xo) =R'^ I UQ dx. (1.7) 



The critical time is defined as in (1.5); the function H{'i]) is defined as Kr]^/^^~'^^ for rj < 1 while 
H{r]) = Krj^'^^ for > 1, with K = K{m,d). Note that for < t < the lower bound means 
'u{t,XQ) > k{m,d){t/ R?Y/^^~^^ which is independent of the initial mass. 

Eliminating the time T . A natural question is to try to recover this sharp results of the good fast 
diffusion range via the present methods. If one wants to do that, one needs upper estimates for the 
minimal life time, that is upper estimates for the extinction time for the MDP, in terms of the L^-norm 
on the ball Br^. We prove the following result. 

Theorem 1.3 Let rric < rn <\. Then, (i) We have sharp upper and lower estimates for the extinction 
time for the Dirichlet problem on any hall Br of the form: 

^ IL, 7?2-d(l-m) ^ rp ^ II ||l-m T}2~d{\~rn) /-i o\ 
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(ii) In that range of m the lower estimates of Theorem imply the lower Harnack inequalities of 
[111201 [211 [H], in the form 



u{t, Xq) > C 



'm,d 



(1.9) 



for any < t < t* and any x G Br, where t^ is given by (1.26). 



This result shows that the form of the lower bounds given in Theorem |1.1| is sharp, since it allows to 
obtain sharp local lower bounds not only in the good fast diffusion range. And it also applies in the 
very fast diffusion range, that is the new interesting part of this paper. We are thus led to the question 
of eliminating all extinction times from the estimate, i.e., replacing T or Tm by some information on 
the initial data, also in the range < m < nic- 



Theorem 1.4 Let < m < rric and let u he the solution to the FDE under the above assumption that 
uq G L^^^{M.'^). Let Xq be a point in O and let d{xo, 90) > 3R. Then, the following inequality holds for 
alio <t <T 



with Cl and C3 given positive constants depending on d. 



u"'{t,xo). 



(1.10) 



We can also obtain formulas in terms of the norms ||tto||LP(_Bfl(xo) P ^ ^'c, that can be seen 

below. 



Obstruction to a simpler estimate with norm 

The presence of the extinction time T in the lower estimates, or equivalently of some norm of the 
initial data, is a drawback in the formulas that is not present in the original Aronson-Caffarelli estimate 
for m > 1, or in the version of the authors for m € {rric, 1) in the whole space. But it is a consequence 
of the 'bad' behaviour of the fast diffusion equation for low values of m, a fact that can be seen in 
different ways. 

Thus, we will show here that the local lower estimates cannot depend only on the local norm of 
the initial data when < m < nic- We do it by means of a counterexample based on the behaviour 
of solutions with data that approximate a Dirac delta. We solve the FDE for smooth and positive 
initial data ^p{x) E L^(M'^) with integral equal to 1. We assume that ip is radially symmetric, compactly 
supported and decreasing with |x|. We obtain a smooth and positive solution u{t,x) defined in a 
cylinder Qt^ and vanishing identically at some t = Ti. The scale invariance of the equation implies 
that the solution corresponding to data (pk{x) = k'^(p{kx) is 

Uk{x) = k'^u{k'''t,kx), a = d{l-m) - 2>0, (1.11) 

so that it has extinction time T^ = Ti k" . As /c ^ oo it is clear that Mfc(0,t) converges to the Dirac 
delta. We also observe that Tk — > cxd, so that we lose the previous estimates. On the other hand, we 
see that losing the estimates is inevitable. If we consider a point xq very close to x = and take a 
radius R > \xq\, then ||tifc(0, a;)||Li(^^(^,Q')) = 1. However, by continuity of u with respect to the initial 
data at t = 0, X large, we have 

Uk{t,xo) = k'^u{k~''t,kxo) 
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(note that Jj^^ Uk{ti x) dx < 1 at all times). This means that no lower estimate could be uniformly valid 
for this sequence. 

A scaling argument was used by Brezis and Friedman [11^ to prove that there exist no weak solutions 
with initial data a Dirac delta. 



1.1 Positivity for a "minimal" Dirichlet Problem 

We will assume that < m < 1 in the study of positivity (cf. the comment in the Introduction). Since 



m > we eliminate the factor 1/m from equation (0.1) for simplicity without loss of generality. As a 



preliminary step, we first prove positivity for a problem posed on a ball of radius i?o, zero boundary 
data and particular initial data. Since the problem of getting quantitative positivity estimates has been 
successfully studied in |8J in the range nic < m < 1, the techniques we introduce are mainly aimed at 
producing positivity in the cases < m < rric, where previous methods failed. 

Specifically, we shall consider the following Dirichlet problem on the ball S/Jq C M'^: 

' dtu = A(7x™) in = (0, T) X 

u{0,x)=uq{x) in Bfj^, and supp(no) ^ i?/? (1-12) 

u{t, x) = for t > and x £ dBfj^ , 

where Rq > 2R > 0. We only consider nonnegative data and solutions. The problem admits a unique 
solution u £ C([0,cxd) : L^(Bjiq)) for every uq £ L^(^Bji^^, We will refer to this problem as the 
minimal Dirichlet problem, or more briefly, the minimal problem, because obtaining positivity for 
solutions to this problem implies in an easy way local positivity for any other problem, thanks to the 
comparison principle. The solution vanishes in finite time; let T > be the finite extinction time, 
shortly FET. Later on we would like to eliminate the dependence of the results on T and make the 



estimates depend only on the initial data, see Section 1.4 

Our goal is to obtain positivity with a quantitative estimate for this "minimal" problem. Our most 
novel idea consists in passing the information on the initial data via the flux of the solution on the 
boundary of the ball B2R into an averaged positivity result outside the ball for times that are not 
too small, more precisely on the annulus Aq := Br^ \ B2R. This property can be interpreted as the 
expansion of positivity outside a ball in which the initial datum has nonzero mean. It is in some 
sense it is analogous to the expansion of positivity already introduced by DiBenedetto et al. , see e.g. 
|19| ITSl I2T] for the upper m-range . The expansion of positivity turns out to be a key tool in proving 
lower Harnack also in our case. 

Once we have proved that positivity spreads out from a ball, then for suitable positive times the mean 
value of the solution on an annulus is positive. We then "fill the hole" in the middle using Aleksandrov's 
Reflection Principle, cf. the Appendix and [8]. In this way we arrive at the positivity result in the 
inner ball for any positive time. 



1.1.1 Flux and transfer of positivity 

We start the proof of the positivity results for the minimal problem by a result on mass transfer 
to an outside annulus based on the flux across an internal boundary. We recall that Rq > 2R and 
^0 '■= Bfi^{xo) \ B2r{xq). In order to simplify the final formulas, we write A = Ro/2R > 1 (we take for 
instance Rq = 3i2). 



8 



Lemma 1.5 (Flux Lemma) If u is a positive smooth solution of the Minimal Problem (1.12) in Qt 
with extinction time T > 0. Then, the following estimate holds true 



ko{Ro-2Rf f u{s,x)dx< f [ u'^dxdt, (1.13) 
Jbr^^ Js Jao 

for any < s < T, and any < 2R < Rq, and for a suitable constant = ko^d). 

Proof. We shall use a C°° test function ip{x) that is supported in the ball Bji^ and takes the value 
1 in B2R- It is clear that we can choose ip such that there exist a constant kQ > depending only on d 
such that 

l^"'^)! ^ (Ay • 

Let < s < t < T. We compute 

/ / dtu^dxdt= j j A('u'")(/?dxdt = / / u™A(/?dxdi+ / / du{u'^)^dadt 

J s J Ao J s J Aq j s j Aq j s j dB2R 

+ f f [du{u'^)^-u'^du^]dadt- f I vJ^dy^dadt. 

Js JdBRg Js JdB2R 

We remark that the last three integrals vanish since 99 and u = vanishes identically near the boundary 
dBpiQ, and d^ip = on dB2R. We also have 

/ dtUifdxdt= I u{t,-)ipdx — / u{s,-)(pdx 

S J Aq j Aq j Aq 

Hence, 

u{t,-)ipdx- I u{s,-)(pdx= I I u^l\ipdxdt+ j f d,y{u"^)ipda dt. 

Aq j Ao Js j Aq Js J dB2R 

We will use this equality with t = T, T = T(uq) being the finite extinction time for the solution to 



Problem (1.12), so that we obtain 



u{s,-)Lpdx = - I I u"'Aipdxdt+ [ I du*{u'^)^dadt, (1.15) 

Ao Js J Ao Js JdB2R 

where v* is the exterior normal to B2R., which is the opposite of v which is the exterior normal to the 
inner boundary of so that du*{vJ^) = —djy{u^). 

On the other hand, a simple calculation shows that 

/ {u{t,x)-u{s,x))dx= f j dtudxdt= j j A(u")dxdt=/" / d^,*{vJ^)da dt. 

J B2R J s J B2R J s J B2R J s J dB2R 

Letting t = T, with T as above, we obtain 

u{s,x)dx= I I d^*{u"')dadt. (1.16) 

B2R Js J dB2R 
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Joining equalities (1.15) and (1.16) we get 



u{s,x)dx= / u{s,x)dx+ / u{s,x)dx 
Brq J B2R J Ao 



s J Ao 



We conclude by using estimates (1.14) for Ay?: for any < 2R < Rq we then get 



u(s, x) dx < 



u(s, x) dx 



Rq 



[ L '''' ^ [ L 



The proof is complete. □ 



Remark. Lower Bound on the Extinction Time. As a first consequence of this Lemma we can easily 
obtain useful lower estimates for the FET: 



ko{Ro-2Rf [ u{s,x)dx<f [ dx dt < (T - s) Vol(^o) / u"'{s,x) 

Jbr^ Js Jao J Ao 

m 



dx 



<{T-s) Vol(^o) 



m(s, x) 



dx 



Vol(^o 



Vol(^o) 



<{T-s) Vol(.4o 



1 1— m 



u{s, x) dx 



B 



Ro 



where in the first step we have used the mean value theorem for the time integral (see details in Step 
2 of next section), with s G (s,r), in the second step the Holder inequality, and in the third step we 
used the contractivity of the global L^(Sijg)-norm. Letting then s = gives the desired lower bound, 
once we notice that supp(tio) ^ Br 



ko {Ro - 2Rf 



I Br ^0 



Vol(Ao) 



1— m 



< T. 



^1.18) 



1.1.2 Pointwise lower estimate for initial times 

We have just shown how positivity of the initial datum propagates on the annulus in the weak form 
of a positive space-time mean value. We will now see that this is sufficient to fill the hole inside the 
annulus. As in the study of the exponent range rric < m < \ performed in [S], the estimate uses a 
critical time that is defined in terms of the initial norms. In the present case it is given by 



{Ro - 2Rf 



Ibr ^0 dx 



Vol(Ao 



l—m 



(1.19) 



where ko as in the Flux Lemma 1.5 Note in passing that the positivity result that follows, formula 



(1.25), implies that this quantity is less than T. 



Obtaining the lower bound needs several steps. 
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• Step 1. Time Integrals. Holder's inequality, together with the fact that the global (i?/jp)-norm 
decreases, gives 



Ao 



u{t,x)"'dx < Vol(ylo) 



< 



Vol(ylo) 



1— m 



1— m 



u{t, x) dx 



Ao 



< 



Vol(Ao) 



1— m 



u(t, x) dx 



B 



n(0, x) dx 



Vol(^o) 



1— m 



Uq dx 



Br 



since supp(Mo) ^ Bn . For any < s < t we then have 



s J Aq 



i(r,x)™dxdT < Vol(^o) 



1— m 



uo dx 



Be 



(t-s) 



We use this estimate together with estimate (1.13) to get 



ko{Ro-2Ry uo{x)dx< / n"dxdt 
Jbh Jo Jao 



u"^ dxdt + 



JAo 



u"^ dxdt 



U JAo 



< Vol(^o)' 







m r 


m 


I uq dx 


U+ I / u™ dx dt 




-JBr 


Jt, JAo 



[1.20] 



In view of the definition of t^, we can eliminate one term and get 



k2{Ro-2Rf [ no(x)dx < / / u'^dxdt 
J Br Ju J Ao 



(1.21) 



with /c2 = kQ/2. In particular, this means that the left-hand side remains strictly positive. 
• Step 2. We introduce the function 



Y{t) = / 'u™(t,x)dx, 
JAo 

and apply the mean value theorem -for the time integral- to prove that there exists ti G such 
that jf^ Y{t) dt = (T -Q Y{ti) . In other words. 



/ / u"'{t,x)dxdt= {T -U) I u'"(ti,x)dx. 

Jtt J Ao J Ao 



Using now the estimate obtained in the previous step, we conclude that there exists ti G \^^, , T), such 
that ^ 

k2{Ro-2Rf I uo{x)dx< / u"'{t,x)dxdt = {T -U) n"(ti,x)dx, 
Jbr Ju J Ao J Ao 

and this implies that for some ti G [t* , T] we have 

k2[Ro-2Rf f ^^^^^^^^ k2{Ro-2R)' j- f (122) 



T 



Br 



{T - U) 



Br 



Ao 
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Step 3. Aleksandrov Principle. Positivity at the critical time. We can now use the Aleksandrov 



Principle to deduce positivity at xq from inequality (1.22). In fact, 



Ao 



u'"(ti,x) dx < Vol(Ao)'u™(ti,xo) 



:i.23) 



where xq G M'^ is the center of the ball -Bhoj since we know from Aleksandrov Principle, that u{t, x) < 
u{t,Xo) for any x £ Aq and any t > (see Appendix for details). 

Joining inequality (1.22) and (1.23) we obtained that there exists a ti G , T) such that 

k2 (Ro - 2Rf 



Vol(ylo)T Jb^, 



UQ{x)dx < u"^(ti,xo) 



(1.24) 



• Step 4. Positivity backward in time. The last step consists in obtaining a lower estimate when 
< t < ti. This argument is based on Benilan-Crandall's differential estimate, cf. |3] : 



dtu{t, x) < 



u{t, x) 
'1 — m)t 



that is valid for all nonnegative solutions of this initial and boundary value problem. It easily implies 
that the function: 

1 

u{t,x)t 

is non-increasing in time, thus for any t £ (0, ti] we have that 



u{ti,x) < t'^ tl-^u{t,x) < u{t,x) 



since ti<T . It is now sufficient to apply inequality (1.24) to the l.h.s. in the above inequality to get 
k2 (i?o - 2Rf 



UQ{x)dx< ti'"(ti,xo) < t'^ u'^{t,xo) 



1.25 



voi(Ao)r 

This is the inequality we were looking for. 

• Step 5. In order to simplify the final formulas, it is convenient to make a choice for the ratio 
A = Rq/2R > 1 (for instance i?o = 3R). The formula for t^, becomes 

i,=c'oi?2-<i(l-™)||^^||l-rn^^ (1.26) 

and Cq > depends only on d and A. We have proved the following positivity result. 



Theorem 1.6 Let < m < 1 and let u be the solution to the Minimal Problem ( |1.12 ) and let T 
T{uq) be the MET. Then T > 2t^,, and the following inequality holds true for all t G (0,t=K] 



u"'{t,xo) > c\R'-''\\uohiiBj,)T'^ t-- 
where c'l > depends only on d. 

For the particular time t = we get 



(1.27) 



u{t^,xo) 



§j^^u{0,x)dx ^ 



yc'^iR^Ty/^^-'^'' 6 u{0,x)dx. 

J Br 
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1.1.3 Estimate of Aronson-CafFarelli type for Very Fast Diffusion 



It is interesting to present the above result in the form that has been used by Aronson and Caffarelli 
in their work [1]. We have to argue as follows: we have arrived at the following alternative 

either t* < t or ^ — / uo{x) dx < t i-'"u™(t,xo) 

curf(A^- l)i?'^-2rT^ Jbr 

Writing the expression of t*, we either have 

u^diX'' - 1) 



[ Uoix) dx < Cl t , Cl 



or 

R-" I uo{x)dx<C2TT^R-^t-T^ u'^{t,xo), C2 = 5^4t^^ 

J Br k2 (A^ - 1) 

We now sum up the two expressions to get 



[ uo{x)dx <CiR-'^'^^-"'h^ +C2T^R-^ vJ^it.xo). (1.28) 

JBn 



I Br 

with Cl and C2 given constants depending on d and A > 2. 



1.2 Proof of Theorem 1.1 and Corollary 1.2 



The previous results will now be used to prove uniform positivity on balls for any local solution of the 
FDE problem. We proceed as follows: let u be a positive and continuous weak solution of the FDE 
defined in Q = (0,T) x O taking initial data n(0,x) = uq{x) in VI. We make no assumption on the 
boundary condition (apart from nonnegativity) . Let us select a point xq & Q W^. Select two radii 
i?o ^ 3i? > so that Bji^^{xq) C 17, that is Rq < d\st{xQ,d^) . In the case J7 = M'^ there is obviously 



no restriction on Rq. Let U£)he the solution to the corresponding MDP, as defined in (1.12). It has 
extinction time Tm- By parabolic comparison, it is then clear that U£) < u in Q = (0,T') x i?/jp(xo), 
where T' = min{T, T^}, hence we can easily extend the positivity results for the MDP obtained in the 



previous section to any other local weak solution, either in the form given by Theorem 1.6, or in the 



Aronson- Caffarelli form (1.28). This concludes the proof of Corollary 1.2, and as a consequence we get 
Theorem |1.1| □ 



1.3 Lower estimates independent of the extinction time for rric < m < 1. Proof of 
Theorem [TT3] 

The proof is divided in some short steps. Here, rric < m < 1. 

• Reduction. Let uii{t,x) be the solution to the homogeneous Dirichlet problem on the ball Bn, 
corresponding to the initial datum uq G 'L}{Br) and having extinction time T{R,uo). The rescaled 
solution 



u{t, x) 



M 



t X — Xq 



j^2-d{l-m) j^l-m 



M 



uo{x) dx 



(1.29) 



Bf 
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allows us to pass from a solution with mass M defined in the ball of radius R centered at xq to a 
solution u with mass 1 in the ball of radius 1 centered at 0. The extinction times have to be rescaled 
accordingly, T{u) = R^-d{i-m) j^i-rmf ^ ^j^^^^ if ^ j,^^^ Therefore, we will work with rescaled 
problems and solutions. 

• Barenblatt solutions. We now consider the solution B of the Dirichlet problem posed on Bi, 
and corresponding to the Dirac mass as initial data, i3(0, •) = (5o, and zero boundary data, that we call 
Barenblatt solution. First we recall that by approximation with L^-data, or by comparison with the 
solutions of the Cauchy problem, it is easy to check that the smoothing effect applies and 

Bit, x) < Cm,d t~'^^^ , for any (t, x) G [0, +oo) x Bi . (1.30) 

Moreover, it is known that this is the solution that extinguishes at the later time among all nonnegative 
solutions with the same mass of the initial data and same boundary data. Such comparison is a 
consequence of the concentration-comparison and symmetrization arguments developed in detail in 
|32| |3T] . Thus, we need to prove that the Barenblatt solution extinguishes in finite time T. The proof 
is based on the fact that it is bounded for t > to > 0. 

• Solution by separation of variables. Consider now the solniion Ur{t,x) = S{x){Ti — t)i/(i-™') 
of the Dirichlet problem on 5^, r > 1. It corresponds to the initial datum Wr(0,x) = S{x)tI^^^'^ ™\ 
and extinguishes at a time Ti, to be chosen later. Here, S is the solution to the stationary elliptic 
Dirichlet problem AS"" + 1/(1 — m)S = on -Bij^, and therefore it can be chosen radially symmetric, 
S{x) = S{\x\). It will also be nonincreasing in r = By standard regularity theory S'(x)'" can be 
bounded from above an below by the distance to the boundary. The parameter Ti, extinction time of 
lAr can be chosen at will. To fix it, we pick a time to > and define the Ti through the relation 

S(l)(ri - to)^ = Wto'''''- (1-31) 



• Comparing the two solutions. We now consider the homogeneous Dirichlet problem on [tQ,T\ x 
Br, and we compare the Barenblatt solution B with the solution Ur constructed above in the cylinder 
Q = -Bi(O) X [to,Ti). At the initial time to we know by construction that B{tQ,x) < Ur{tQ,x). The 
comparison of the boundary data is immediate. By well-known parabolic comparison results, this 
implies that B{t,x) < Urit,x), on [toiT^i] x Br, and hence the extinction times satisfy 



T <Ti 



5(1) c 



1— m 



+ to. 



(1.32) 



We only need to choose a to G (0, T) to obtain an expression for the upper bound of T that depends 
only on m and d (the reader may choose to optimize the expression for Ti with respect to to). 

• Conclusion. As a consequence of the above upper bound, we know that any solution to the Dirichlet 
problem on the unitary ball and with unitary initial mass extinguish at a time T < T < T{m,d). 
Rescaling back to the original variables we have proved that any solution u/j to the Dirichlet problem 
on the ball Bpi and with initial mass M = J^^ uq dx extinguish at a time Tr that can be bounded 
above by 



The lower bounds come from the fact that t^, <T and is given by (1.26). 
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• The lower Harnack inequality. Inequality (1.9) follows by plugging the upper bound (1.8) 
into the lower bound (1.27). □ 



The reader should notice that the properties that we have used are typical of the good fast diffusion 
range, rric < m < 1, and cannot be extended to the very fast diffusion range, m < rUc- 



1.4 Lower estimates independent of the extinction for < m < rric. 

The presence of the minimal extinction time Tm = T in the formula for the lower Harnack inequality 
responds to an essential characteristic of the problem. Actually, lower estimates in terms of only 
norms cannot be true for m < rric as we have shown at the beginning of this section: there is no 
positive lower bound at a time to > and a point xq that depends only on tQ,R and the mass of uq 
in Br{xq). Similar examples can be constructed if uq G L]'^^(M'^) with p < pc, and we can be found in 
[3T] . Chapters 5 and 7. 

Fortunately, controlling the local (or global) norm gives a control on the MET T, and in this way 
we get valid lower estimate without T, as we explain next. 



1.4.1 Estimates in terms of the L^^ norm 

In [3] Benilan and Crandall prove that for any < s < t < T, and for any m < rric 

Mml:"^ < MsniT"" - jc^^ (t-s), with /c^, 

where 52 in the constant of the Sobolev inequality 

|2- < 52IIV/II2 



8 [d{l -m)-2] Si 
(d - 2)2(1 -m) 



:i.33) 



(1.34) 

and the above estimate holds for any solution with initial datum uq E U''^. We also stress on the fact 
that the constant Kp^ is universal in the sense that it only depends on m and d. As a consequence of 



(1.33) , we have the following universal upper bound for the extinction time 

T{uo) < JCp^ \\uo\\ll 



\l—m 



(1.35) 



We remark that while for lower bounds on FET we only need local information on the initial datum, 
upper estimates for the FET require global information. Fortunately, in the minimal problem that we 
are considering, global and local are equivalent since uq{x) = for |x — xo| > R- 

Proof. We sketch here the proof for the reader's convenience. It is well known that the time derivative 
of the global norm of the solution u{t) of the MDP problem under consideration is given by 



— \ 



n{t)\\l 



< 



4p(p — 1) 
{p + m — 1)2 

{p + m — 1)2 



p-\-m — l 

Vu 2 



dx 



(p+m-l)2' 

u 2 dx 



4p{p- 1)5| ..p+m-l 

■(p + m-l)2 ' 



^1.36) 



where in the last step we used the Sobolev inequality (1.34) applied to the function / = u^p^^ ^)/'^^ 
where 2* = 2d/{d — 2) and 52 is the Sobolev constant. 
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Notice that if m > rric, then pc < 1, so that the global L^-norm increases, and this originates the so 
called Backward Effect, see e.g. [31]. This explains our assumption m < me - Moreover 



Pc + m-l=pc[l 



{pc + m- 1)2* 



Ap,{p,-l)Sl ^ 8 [d{l - m) - 2] Si 
[Pc + m-iy (d - 2)2(1 -m) 



> 0, 



so that (1.36) becomes 



d 



8[d(i-^)-2]5|„ p. 



^^u<t)¥p:< (d_ 2)2(1 _m) 



integrating the differential inequality gives the bound (1.36) for any < s < t. Letting s = and 



t = T{uo) in (1.36) finally gives (1.35). □ 



Application to Theorem 1.4, Using this bound, we can now formulate the second version of the 
lower Harnack estimate, reflected in the theorem. The proof is immediate when uq is as in the minimal 



problem, since in that case local and global norm is the same. Comparison as done in Subsection 1.2 



allows to pass to the general solutions. Notice that in this way we use a local L^" norm, not the global 
one! 

Remark. When we have not only the Sobolev inequality, but also the Poincare, we can prove similar 
estimates for any m E (0, 1) . This happens for instance for problems posed in bounded domains, or 
for the minimal Dirichlet problem. 

1.4.2 Estimates in terms of other norms 

Proposition 1.7 Let m<l,a>l,R>0 and let u be the solution to the Dirichlet problem 

ut = ^A{u"') m{0,T)xBaR 
u{0,x) = uo{x) in BaR, and supp(tio) ^ Br 
u{t, x) = for t > and x £ dB^R 

with Uq G L^^BaR), with p > max{pc, 1} = max{(i(l — m)/2, 1} . Then the following estimate 

I l—m II / \ 1 1 1— m 



\u{t)\\ 



|^^(s)||p < -/Cp (t - s) 



(1.37) 



hold for any < s < t , where 



JCp = 4 



2pc 



[p + m — l] 

and where S2 is the Sobolev constant ofW^ and V is the Poincare constant on the unit ball. 

1 2 / \ -* 

Proof. First consider, for any / G Wq' [BaR) , the Sobolev and Poincare inequalities: 

I2. < S2 IIV/II2, and II/II2 < VaRWVfh , 
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where 2* = Idjid — 2) , and where the constants 52, the optimal Sobolev constant on M'', and 7-", the 
Poincare constant on the unit ball, only depend on the dimension d . By combining them through the 
Holder inequality, we then get for any q G (2,2*) 



|/|L<I|/Il2~l/Ilt <[^«^]'~''5|||V/||,. 



Let now 



p+m — 1 



2p 



and 



d(l — m) Pc 



p + m — 1 2p "p 

We remark that g < 2* if and only if p > pc, while q> 1\i and only \i q < oo . We obtain then 



\u 



2pc 



p+m— 1 


2 


p+m — 1 


Vn 2 


:= /Co 

2 


Vu 2 



The derivative of the global L^-norm then satisfies 



dt 



4p(p — 1) 
[p + m — l] 



p+m — 1 

Vu 2 



2 4p(p-l)/Co-^ || ||Ji-i 
< ^ m L 



[p + m — 1 



2 II hp 



(1.38) 



(1.39) 



where in the last step we used (1.38). Integrating the differential inequality over [s,t] C [0,T], gives 



|l-m II ,M|l-m_^ 4(1 - m)(p - 1) ^,_2('l-2£^ 



\u{s)\\ 



[p + m — l] ' 



[V aR] 



W ^2 Mi - ^) • □ 



Upper Bounds on the Extinction Time. The above estimates (1.37), prove that any solution of 
the Dirichlet problem extinguish in finite time, and this is not surprising, but they also provide an 
Upper Bound for the extinction time T , indeed letting s = and t = T, we obtain 



T < /C„ ^lluo 



1 1—m 
\p 



2pc 



4(1 - m)(p - 1) ^ ^ ^ " "P 



Notice that in the limit p ^ pc we recover the previous result (1.33). Summing up, the above result 
proves that a global Sobolev and Poincare inequality provides that the solution extinguish in finite time 
T and an gives a quantitative upper bound for T . 

Remarks. These results can be extended to different domains or manifolds in a straightforward way, 
the only important thing is to have global Sobolev and Poincare inequalities, as already studied by the 
authors in [7^ , in the case of Riemannian manifolds with nonpositive curvature 

Using this bound, we can now formulate a version of the lower Harnack estimate similar to Theorem 



1.4 We leave the easy details to the reader. 



2 Part II . Local upper bounds 



In the second part of this work we turn our attention to the question of upper estimates for solutions 
with data in some Lj'^^, p > 1, and obtain quantitative forms of the bounds that are sharp in various 
respects. The range of application is all m < 1, even m < 0. We assume moreover that d >3, which is 
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the interesting case also for the lower estimates, in order to avoid technical complications which break 
the flow of the proofs and results, but we remark that the qualitative fact, the existence of local upper 
bounds, is also true for d = 1,2. 



As a preliminary for the main result, we devote Section 2.1 to establish the conservation of the local 



integrability of the solutions and the control of the evolution of the local norm for suitable p > 1. 
Let u = u{t, x) be a nonnegative weak solution of the FDE for m < 1 defined in a space-time cylinder 
Q = (0, T] X for some Ro, T > 0. This is the form of the estimate we get: 



Br{xo) 



\u{t,x)\Pdx 



(1— m) /p 



< 



\u{s, x)\P dx 



(l-m)/p 



+ K{t-s), 



for any Rq > R and < s < t < T. It is valid for all m < 1 if p > 1, p > 1 — m. The dependence 
of X on i? and Rq is explicitly given in Theorem 2.3 below. The estimate extends Herrero-Pierre's 
well-known estimate to p > 1 and is valid for m < 0. 

The main result of this part is the local upper bound that applies for the same type of solution and 
initial data, under different restrictions on p. Here is the precise formulation. 



Theorem 2.1 Let p > 1 if m > rric or p > pc if m < rric- Let u he a local weak solution to the FDE 
in the cylinder (0,r) x i7 C (0,-|-oo) x W^. Then there are positive constants Ci, C2 such that we have 



uit,xo) < 



|no(a;)|^ dx 



+C2 



where Rq < dist(xo,9$^) and the constants Ci depend on m,d and p. 



Rl 



(2.1) 



We recall that t?p = l/(2p — d{l — m)) = l/2{p — pc). Note that the constants Ci do not depend on 
the radii, but only on m, d and p. An explicit formula for them is given at the end of the proof, but 
we point out that such values need not be optimal. The result is proved in Sections 2.2 and 2.3, A 
similar smoothing effect result has been proved for the first time by Herrero and Pierre in [2^ in the 
good fast diffusion range rric < rn <\ using p = 1, but it is new in the range m < where HP's result 



cannot hold in view of solutions like (0.2 ). HP's technique relies on stronger differential estimates that 



do not hold in the subcritical fast diffusion case or on the local setting; our impression is that their 
techniques can not be adapted to the very fast diffusion range. Related estimates for p> 1 are due to 
DiBenedetto and Kwong, [19^ , and Daskalopoulos and Kenig, [14], but as far as we know no results 
cover the very fast diffusion range. Finally, note that the smoothing effect L^^^ into is false for 
exponents p < pc &s, has been demonstrated in [31j. In fact, that monograph studies the existence of 
the so-called backward smoothing effects that go from L*'(M'^) into L^(M'^) for p < pc- 



The local bound in (2.1) is expressed as the sum of two independent terms, one due to initial data. 



the other one due to effects near the boundary. The estimate is optimal in the following senses: 

(i) The first term responds to the influence of the initial data and has the exact form that has been 
demonstrated to be exact for solutions that are defined in the whole space and have initial data in 
LP(M'^), see |31| Chapters 3,5]. By exact we mean that the integral is the same (but extended to the 
whole space) and the exponents are the same, only the constant Ci may differ. We can then recover the 
global smoothing effect on M*^, just by letting i? — > c« so that the second term disappears; as mentioned 
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above the constant Ci is not the optimal one: the best constant for the global smoothing effect on R*^ 
has been calculated by one of the authors in |31j . 

(ii) The last term accounts for the influence of the boundary data and is special to fast diffusion in the 
sense that it does include any information on the precise boundary data, thus allowing for the so-called 
large solutions that take on the value u = +00 on OBr. The term has the exact form prescribed 



by the explicit singular solutions (0.2). This last term has the meaning of an absolute bound for all 
solutions with zero or bounded initial data; thus it can also be interpreted as a universal bound for 
the influence of any boundary effects. In applications it is interpreted as an absolute damping of all 
external influences. 



2.1 Evolution of Local L^-norms 

A basic question in the existence theory is obtaining a priori estimates of the solutions in terms of the 
data measured in some appropriate norm. The peculiar feature of the FDE is the local nature of the 
estimates. A fundamental result in this direction is the local Lioc'Lioj. estimate due to Herrero and 
Pierre (which is valid for m > 0): 

Lemma 2.2 Let u,v ^ C ([0,-|-oo) ; Lj'Q^(]R'^)) be weak solutions of 

dtu = A{u'^/m), < m < 1. 

Let R > 0, Rq = XR with A > 1, and xq G M'^ be such that Bi = (xo) C M°'. Let moreover v < u 
a.e. Then, the following inequalities hold true: 



Br. 

for any t,s >0, where 



[u{t, x) — v{t, x)] dx 





1— m 




dx 


< 


[/. 









H 1—m 



Cl 



{Ro - R) 

and the constant ci > depends only on m, d. 



[u{s, x) — v{s, x)] dx 



+ Kr^r„i \t - s\, (2.2) 



(2.3) 



This result was proven in Prop. 3.1 of [23] and has been generalized to the case of fast diffusion on a 
Riemannian manifold by the authors in [7] . Our goal here is to extend such a result into and Lj^c'-'^ioc 
estimate for suitable p > 1. This estimate has two merits: first, it is valid for all —00 < m < 1; second, 
it is needed for some values p > Pc for the proof of boundedness estimates. 



Theorem 2.3 Let u G C ((0, T) ; Lj'Q^(O)) be a nonnegative weak solution of 

dtu = Aiu'^/m), 



(2.4) 



and assume that u{t, ■) £ Lj^^(i7) for some p > 1, p > 1 — m, and for all < t < T. Here, 0, is a 
domain in M'^ that contains the ball Bi = i?R„(xo). Then, the following inequality holds true: 



\u(t, x)|^ dx 



(1— m) /p 



< 



|n(s, x)|^ dx 



(1—m) /p 



+ Kr^Ro^p {t- s), 



(2.5) 
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for any < s < t < T, where 



PCm,d 



and the constant Cm,d > depends only on m, d. 



Yol {Br, \ Br) 



(1— m) /p 



> 0, 



(2.6) 



Remarks, (i) The result implies for those values of p that whenever u{s, •) G L|^^(r2) for some s > 
p 

-'loc 



then u{t, •) G L^q^(O) for all t > s. Note that the dependence of the local norm is again expressed 



as the sum of two independent terms, one due to the initial data, the other one due to effects near the 
boundary. 

(ii) Note that the times t and s must be ordered in this result, a condition that is not required in 
Lemma 12. 2[ 

(iii) The last term may blow up as we approach the boundary of 0, (where no information on the data 
is used). Indeed, the constant can be written in the form 



Kr,Ro,p = pCmXo" ""'^'"FiR/Ro), F{s) 



(1 _ gd-^(l'ni)/p 



Now, if xq G we may take Rq = d^XQ, dO,) and R = Rq{1 — e). In that case the constant in the last 
term behaves as e — > in the form 

Kr,r,,p ~ i?o^P~'''^^''e"^' /? = 2 - (1 - m)/p = {2p + m - I) /p. 

(iv) The constant blows up in the limit m ^ 1~ , and this is perfectly coherent, since a similar estimate 
is false for the Heat Equation. 

Moreover, the constant i^_R,i?Q,p, blows up when p oo, thus it does not provide local stability, 
while it provides local stability, for p > pc- 



Proof of Theorem 2.3, (i) Let n > and take a test function -0 G C^(^) and ^ > 0. We can 
compute 



dt 



ipu^ dx = p ijjuP ^dtu dx = —p / \/ {4)vP ^) • V 



m 



dx 



W • (^x^+^-^Vn) dx + (p - 1) / V'"^"^""^|V'u|2dx 



1 



p + m — 1 



VV-- V(uP+™-i)dx + 



4(p-l) 



V 



< 



p + m — 1 
p 



(A0) dx 



[p + m — 1)^ 
Ap {p — 1) 



{p + m — 1) 



[ ^ 


p+m— 1 

Vn 2 


2 

dx 


Jn 






p+m — 


1 2 





(2.7) 



dx 



p + m 



- [ I A^/; I mP+"^-Mx. 
1 Jn 



This computation holds true for any p > 1, and any m G M, when one replaces, in the limit m — > 0, the 
quantity {u^)/m with logu; we also have to replace u^~^"^~^ / {p + m — 1) by log('u) if p + m — 1 = 0. 
Of course, when p + m — 1 < the last term may be infinite, since it contains -^^+'"-1 so we make the 
assumption p > 1 — m. 
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(ii) Under such assumptions, inequality (2.7) implies that for any solution u>0 we have 



/ i;{x)u{t,x)Pdx < ^— y / |A(^(x))| 



< 



|p + m — 1| 



1 — m 
P 



ip{x) u{t, xY dx 



p+m — 1 



'i/'(x) x)"^ dx 



(2.8) 



where in the second step we have used Holder's inequality with conjugate exponents p/(l — m) and 
p/{p + rn — 1), and where 



p 



p + m — 1 



|A(V'(x))| (i-™)-;/;^ (1-™) dx 



(2.9) 



We will check below that this quantity can be made finite by a proper choice of ip. Assuming this for 



the moment, formula (2.8 1 can be expressed as a differential inequality of the form 

y'(r) < Cy'-%T) 

where ?/(r) = ip{x)\u{T, x)\p dx, C = C{ip) and e = {1 — m)/p £ (0, 1). Integrating such differential 
inequality over (s, t) lead to 

y%t)-y%s)<Ce{t-s) 



that is 



i;{x)\u{t,x)\Pdx 



(1— m) /p 



< 



■0(x)|M(s,x)|^dx 



P 



for any < s < t. This will immediately imply the statement, once we prove the bounds 

(1 — m) 



P 



Ciip) = Krxp < +0O. 



(2.10) 



(iii) We only have to verify the form of the last bound. To this end we consider a function ij) = (p^ £ 
C^{M), with 

0<ip<l, ^p = I in Br, ^p = outside Bxr (2.11) 

with A = Rq/R > 1. Moreover, we will assume that Lp is radially symmetric and ip{x) = ^(|x|/-R) , 
where ^ : M — > M is a C^(M) function such that: 

< ^(s) < 1 , p{s) = l,forO<s<l, ^ = 0,fors>A 

where A > 1 and |x| is the distance from a fixed point. We then have 



p p p 

|A (^(x))| i-'"'i/;(x) = (^(x) L 



b{b-l)ip^-^ |V^|Vfe/~^ Av9 

{b-2)p 



P 

1 — m 



P I, 1 P I \^ — ^}P o 



P 
1 — m 



(2.12) 



21 



the last inequality follows from the fact that we are considering a radial function < '^{x) = (p{\x\/R) < 
1, with b > Then, we compute 

\Vip{x)f = R~^\^'{\x\/R)\^\V\x\\'^ < R-^\ip'{\x\/R)\'^ <c'^ R-'^ 
\A^{x)\ = \R-y{\x\/R)\V{\x\)\'^ + R-^Tp'{\x\/R)A\x\\ 



< 



R 



l'''<l^-|/'''l + |^(H/ii)|!^ 



R 



< 



i?2 ' 



where in the last step we used the fact that Aip is supported in Aji x = Bxr \ and that the smooth 
function Tp has bounded derivatives in Ar^\ 



x|/ii)|<^ = c'. 



^\\x\/R)\ + \Tp"{\x\/R)\< 



-0 _ // 



(A -1)2 

we just remark that this last estimate depend on an explicit choice of the test function Jp. 
Inequality (2.12) together with (2.13) gives 



(2.13) 



' c'i + {d - l)c'^ 



P 1 I {b~2)p 



V 
1 — m 



< Hb - 1)] 



P 

1 — m 



[{\-l)R\ 



P 

1 — m 



[(A-l)i?]i 



2p_ 



if 6 > c'p = 6(6 — 1){c'q + {d — l)co). An integration over Af(^\ gives: 



(1 — m) 



P 



< 



1—m 
p + m — 1 

1—m 



|A(?/;(x))| (1-™) V'^ dx 



/ 



1 — m 
P 



p + m-1 [(A - l)i?' 



[{\-l)R] 



^\o\{Ar^x) ^ ■■= Kr^x,p < +00 



where Cp = 6(6 — l)(co^ + (d— l)co)(l — m)/(p+m — 1), and 6 > 2p/(l — m), we can choose 6 = 3p/{l — m) 
to get 

where ^ is independent of p. The proof is thus complete. □ 



2.2 Smoothing effect in terms of space-time integrals 

In this section we are going to prove a first version of the Local Smoothing Effect for the FDE. More 
precisely, we prove that Lj'^^ regularity in space-time implies LJ^^ estimates, even when m < mc, on 
the condition that p must be large enough. The estimates are local, both in space and in time, but 
uniform on balls and the dependence is quantitative. We consider a nonnegative weak solution of the 
FDE for m < 1 defined in a space-time cylinder Q = (0, T] x Br for some R,T > 0. 

Throughout this section T will not denote extinction time. 
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Theorem 2.4 Let u and m he as above, and let p > 1 if m > rric or p > pc if m < nic. For any two 
finite cylinders Qi C Qo, Qi = {Ti,T] x Bji., with < i?i < Rq, and <Tq <Ti <T, we have 



sup|u| < Cloc 
Qi 



1 



{Ro — Ri 



+ 



1 



d+2 
2p+d(m-l) 



uPdxdt + \ol{Qo) 



2p + d(m-l) 



(2.14) 



Moreover, the constant C\oc depends only onm,d,p. 



The proof presented here uses Moser's iteration process, and borrows some ideas of [13] and [21]. We 
will consider nested space-time cylinders, in order to obtain the first estimates needed to prove local 
SE. The proof will consist of the combination of several partial results, which maybe of independent 
interest, and will be split into several steps. Note that by scaling the proof of this kind of result need 
only to be done for a unit cylinder Qo where Ro = 1 and Ti — Tq = 1, and this is the case that will be 
needed in the sequel. 



Step 1. Space-Time Energy Inequality 

Now we consider a solution u defined in a parabolic cylinder Q = (To,T] x Bji for some R > Ri > 0, 
T > and consider another parabolic cylinder Qi = (Ti,r] x Bft^ , contained in Q, since we also let 
To < Ti < T. Then 



Lemma 2.5 Under these assumptions, for every m < 1 and p > max{l, 1 — m}, we have 



uP{T, x) dx + 





p+m — 1 


// 


Vu 2 


' JQi 





dxdt < C{m,p) 



1 



+ 



1 



{R - i?i)2 Ti - To 



(2.15) 



The result also holds when u is a sub-solution, i.e., ut < Am" 



Proof, (i) We multiply the equation dtu = ^Au"^ by il)"^ u^~^ , with p > 1 to be chosen later in a 
suitable way, we take tp = ij){t, x) any smooth compactly supported test function, and we integrate on 
the cylinder Q = (0, T] x Bji. By definition of local weak solution, we obtain 



w 



4(p - 1) 



p+m— 1 

Vn 2 



-0 dx dt 



{p + m — ly 

We now use Young's inequality: for any "o^, b G M'^, and any 5 > we have 



P+'^-^W-V'Vudxdt. (2.16) 



a ■ b \ < -I a r + 



1 

2^' 



Together with Holder's inequality, this allows to estimate the right-hand side of (2.16): 

4 



iP+"'-'^Vi^-i^Vudxdt 



p + m — 1 



p-'r'm. — l p-\-m — l 

u 2 VV' • ipVu 2 dx dt 



< 



p + m — 1 
2 



1 
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w 



P+^-^lWp dxdt+^ 



p+m — 1 



p-l 



/+-"i|VV|'dxdt + ^^^^ 



p + m — 1 



Vu 2 



p+m— 1 

Vn 2 



ip dx dt 



ip dx dt 
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where in the last step we have chosen S 



p+m- 



j > 0. Putting this calculation into (2.16), we obtain 



iP-^dtuij^dxdt + 



2{p - 1) 
{p + m — 1)^ 



p+m — 1 

Vu 2 



dx dt < 



p — 1 



w 



p+m— 1 



Now, we integrate the first term by parts (with respect to the time variable) 



^dtuxl^'^dxdt 



1 

P JBn Jo 



T 



dt{uP)ij^dxdt 



Br 



uP{T,x)ij^{T,x)dx- [ uP{0,x)^^{0,x) 

JBr 



dx 



Be 



uP{T,x)^'^{T,x)dx - / iiP(0,a;)V'^(0,x)dx 



Br 



1 

P JJq 
2 

P 



vFdt{ij'^)dxdt 
uF ip dt{'4)) dx dt. 



Collecting all the previous calculations, we obtain the first basic inequality: 

2 



1 

P 

+ 



nP(r,a;)V'^(r,x)dx- / ^^(0, a;)V'^(0, 

R. J Br 

2{p-l) 



{p + m- 1)2 JJq 



f)+m — 1 



Vm 2 



dx dt < 



x) dx 
2 



p — 1 



vF ip dtlip) dx dt 
PJJq 

^p+m-l|y^|2 ^^^^ 



(2.17) 



(ii) In order to continue, we assume that the test function ^ satisfies 

• < ip{t, x) < 1, for any (t, x) G Q, and 7/^(0, x) = 0, for any x G Bji 

• ^ = 1 on Qi = [Ti, r] X Br-^ C Q and = outside Q. Of course, we take < Ri < R and 
< To < Ti < r. 



Moreover, on \ we assume that 



IWI < 



R — Ri 



and 



\dti^\ < 



We may then write (2.17) in the form 



p — 1 
P Jbr 



nP(0,x)V'^(r,x)dx + 



2(p-l)' 
(p + m — 1)^ 



< 2 



< 2cj 



1 



p + m— 1 

Vn 2 



IVV'I dxdt + 



+ 



1 



{R-RiY Ti-To 



Ti-To 

Tp'^ dx dt 

p — 1 
P 



u^ip |9t(V')| dxdt 

P+m~l ^ ^p^ ^^^^ 



We observe that 





p+m— 1 




p+m — 1 


// 


Vti 2 


^ dx dt < 


Vu 2 


JJqi 









dx dt 



24 



since Qi d Q and = 1 on Qi, so that we finally obtain 



m,p 



<2ci 



uP{0, x) dx + 
1 





p+m — 1 


// 


Vu 2 


J jQi 





dx dt 



+ 



[w 



p+m— 1 



where 



{R-Ri)^ Ti-To_ 

p-1 2{p-lf 



+ uP ) dxdt 



Cm,p = min 



p ' (p + m — 1)2 



As a conclusion, we have obtained (2.15) with precise constant 

C{m,p) = 2clC~]p . 



(2.18) 



(2.19) 



Note that Cm,p depends also on d though we are not indicating it. 

We conclude by noticing that the proof can be repeated for u sub-solution (with the same regularity) , 
that means, when ut < Au™" the above estimate continues to hold. □ 

Improving the constant. We would like to eliminate the dependence of C{m,p) on p in what follows 
since p will vary (in an increasing way). This dependence takes place through Cm,p- Now, for m > it 
is easy to see that Cm,p > (p— l)/p and we have to assume that p > Po > I, so that C{m,p) is bounded 
by an expression that depends only on pQ and d. 

For m < 0, since we have p > 1 — m we get (p— 1) /p > |m|/(l — m). A lower bound for Cm,p needs the 
last term to be bounded above, and this implies that p must be away from 1 — m, so that we assume 
that p > p'q = {1 + a){l — m) for some a > in which case we get 



Cm,p > min 



1 — m 



,2 1 + 



\m\ 



a{l — m) 



:= C{m). 



(2.20) 



In any case we may write C(m) instead of Cm,p if the family of p's fulfills the stated conditions. 



Final result of Step 1. We need to improve Lemma 2.5 in the following way 

Corollary 2.6 Under the running assumptions, for every m < 1 and p > max{l, 1 — m}, then for any 
Tq < Ti < T , < Ri < R we have 



sup 

se(Ti,T) JBb^ 



n^(s, x) dx + 



p-{-m-~-l 



Vu 2 
1 



{R-Ri 



+ 



dx dt 
1 



Ti-T( 



T 



To J Br 



(W 



p+m— 1 



+ u^) dxdt 



(2.21) 



Moreover, if u is a sub-solution, and u>l, we have 



sup 



u^{s, x) dx + 



T 



Ti JBr^ 



p + Tn — 1 

Vu 2 



dxdt < C(m) 



+ 



{R-Riy Ti-T, 



fl 

J To JBf 



dxdt 



(2.22) 



with C{m) as in (2.20). 
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Proof. First we recall a property of the supremum: there exists a to S (Ti, T] such that 

u^{s,x)dx< I u^{to,x)dx 



sup 

se{Ti,T)JBii^ 



B 



We use this observation together to the result of Lemma 2.5 in two different ways: 



(i) We use Lemma 2.5 with the T replaced by to and still keeping < Tq < Ti < to. We get 
1 



sup / u^{s,x)dx< / u^{tQ,x)dx 

s&iTi,T)JBR^ Jbr^ 

< C{m) 



+ 



(R - Rif Ti - To 



to 



To J Br 



In the sequel recall that to < T. 



(ii) Next, we choose the same Ti and we apply Lemma 2.5 to get 

2 





p+m — 1 


fl 







dx dt < C{m) 



1 



+ 



1 



(R-RiY Ti-To 



To JB 



p+m—l 



(2.23) 



+ uP) dxdt 



(2.24) 



Summing up the two inequalities (2.23) and (2.24) gives the desired inequality (2.21) . 



For the last part, we remark that if we apply inequality (2.21 ) to a sub-solution u > 1, then ^ < 



so that we obtain (2.22) , and the proof is thus concluded. □ 



Step 2. Iterative form of the Sobolev Inequality 

The next lemma is just a different form of the usual Sobolev inequality, adapted to our aims. 



Lemma 2.7 Let f G L'^{Q) with V/ € L'^iQ). We then have 



Ti JB 



f^^dxdt < 2Si 



Ti JB 



f + R^ |V/| dxdt 



sup 

se{Ti,T) 



JB, 



/2('^-l)'?(s,x) dx 



for any a G (1,0"*), and for any < Ti < T and R > , where 



2* 



d—2 m, 

2 , 



ifd=l,2 



and q 



cr 



a* -1 



I, ifd>3 



2, if d = 1,2. 

Here, S2 is the constant of the classical Sobolev inequality \\f\\2* < ^2 (||V/||2 + II/II2); with 2* 
2d/{d -2)ford>3 and 2* = 4 for d = 1,2. 



(2.25) 
(2.26) 



Proof. Since the estimate (2.25) is scaling invariant, it is sufficient to prove it for R = T — Ti = 1, and 
we denote hy B = Bi the unit ball of R'^. By Sobolev and Holder inequalities we then get 



B 



^^dx < 


[ f'dx 


2* 






JB 




Jb 



<2Sl 



V/| dx+ / f'dx 
B Jb . 



sup 

se(o,i) 



/2('^"i)''(s,x)dx 



Bi 



Integrating in time over (0, 1) and rescaling back, gives inequality (2.25). □ 
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Step 3. The Iteration 



In this step we use the inequahties of the preceding steps to start the iteration in the Moser style. 
We first define v{t,x) = max {u{t,x),l]. Then we observe that when n is a local weak solution to 
Ut = A n™' , then t; is a local weak sub-solution to Vt = A w™. It is clear that u < v < 1 + u for almost 
any (t, x) e Q . 

Preparation of the iteration step. Letting = i;P+™--i in the modified Sobolev inequahty 
(|2.25|) gives 



p+m—1 



1 p+m — 1 |9\ 




1 f 


2 1 1 dxdt 




sup — ^ / 

_te(Ti,T) Ri Jbr^ 



^(p+m-l)(a-l)g^^ 
(2.27) 



where Qi = (Ti, T] x Br, C Qo = (To, T] x Br. 

Since we have assumed that v > 1, then uP+™~i < yP and we can use the energy inequality (2.22) 
to estimate the two terms of the right hand side of the above inequality (2.27), in terms of the same 
quantity. First we estimate 



,p+m— 1 



p+m — 1 I 2 



dxdt < 



IIqi 



dx dt 



+ Rj C{m] 
<2RlC{m) 



1 



(i?o — Ri 
1 



)2 + Ti - To 



+ 



dxdt 



{Ro-RiY Ti-To 



v'P dxdt 



Qo 



In the last step we use the fact that i?f C{m) 



1 



+ 



1 



> 1 , which is not restrictive. 



We next estimate the sup term, again using the energy inequality (2.22) , but we replace p with 
{p + m — 1){(T — l)q. If the exponent is larger than max{l, 1 — m} we obtain 



te(Ti,T) Ri 



C{m) 
Rf 



+ 



{R-Ri)^ Ti-To 



Qo 



Summing up, we have estimated (2.27) as follows 
lie 



dxdt < 45| rI ' C(m)^+^ 
yP dxdt 



1 



{Rq — Ri 



+ 



Ti-To 



Qo 



^(p+m-l)(<x-l)g ^^^^ 



(2.28) 



2-2 



Finally, we remark that R^ ' = 1, since - = 2, q being defined as in Lemma 



2.7 



The First Iteration Step. We now use ( |2.28[ ) in the following way: first we choose a e (l,o"*), 

p 



where a* is as in Lemma |2.7[ in such a way that 

{p + m — l){a — l)q = p that is a = 1 + 



q{p + m — 1) 
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A straightforward calculation shows that a € (1, c*) if and only if p > pc . This is where the restriction 
on p appears for the first time. 

We are now ready to begin with the first iterative step, by letting 

Po = P = (p + — — l)g , and pi = (po + — V)o = Po ^1 H — ^ + m — 1 . 



We remark that 



Pi > Po 



Po > Vc 



Estimate (2.28) now becomes 



v^^ dx At < 4cS| C{m) 



1+^ 



1 



+ 



(i(l — m) 



1 



0,1 



1 

dxdt' 









/"/" dxdt 







Qo 



(2.29) 



which is the first iterative step. 
The fe-TH Iteration Step. Letting 

1 



Pfc+i =Pfc(^l + -)+m-l, with pk+i > Pk 
we get the iterative inequality 



Pk>Po> Pc 



yP'^+i dx dt 



Qk + l 



1 

pfe+1 



— ^k,k+l 



v^^ dxdt 



Pfc V lj Pfc+l 



(2.30) 

In order to find a convenient value for Ik,k+i we choose a decreasing sequence of radii i?oo ^ — -Rfc+i < 
Rk < Rq such that < Rk — Rk+i = p/A;^, and a sequence of times < To < < Tk+i —>■ Too < T 
such that Tk+i — Tk = r/k^. This means taking 



p = Ci{Ro - Roo) , T = C2{Too-To 

with ci = 1/ (Efc^o k-^) > 0, and C2 = 1/ fc-^) > 0. Then, 

1+^ 



(2.311 



1 



+ 



1 



{Rk — Rk+iY Tk+i — Tk 



< 45| C(m)^+9 (2 + ^-1) 



(2.32) 



We now calculate the exponents pk- 



Pk+i =Pfc(l + -)+m-l 
1 



1 + 



fc+i 



1 

1 + - 

q 

k+l 

po + - 1) ^ 



fc+i 



1 + 



PQ + {m-l)^ 
I 



n=0 



1 + 



[po - 9(1 



1 + 



k+l 



(2.33) 



+ q{l - m) . 
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notice that 



lim 



1 + ^ 



fc+i 



1 



fc^oo pk+i po + q{m - 1) 

The iterative step now reads 



and 



yP^+i dx dt 



Qk + l 



1 

pfe+1 \l- , 



Now we use (2.32) to estimate 



1 

jPk+l jl' ' 1] Pk+1 



j-l ii Pk+i 
^0,1 



< 



1+ 



hm 



rL 9 J Pfe+i 

-'0,1 



4^ 



1 + 



Po + g(m - 1) 



dx dt 



(2.34) 



X A; "fe+i {k - 1) ffc+i (A; - 2) 



4 1 



2 Pfe+i 1 (2.35) 



i+i 
Jo Jl ' 



Pk+ 



Moreover, passing to the Hmit in (2.34) when A; — > 00 , we get (we refer to Appendix A3 for further 
details) 



sup \v\ < J^o+'Xm^^) jPo+«(— 1) s^e^iQ+i) 



yP" dxdt 



Qo 



PQ + q{m-l) 



(2.36) 



We have estimated the constants to ensure that they remain bounded in the hmit k +00, see the 
Appendix for the details. Moreover, these constants blow up when Roo Rq or Too Tq: indeed, 
while Jo and si only depend on m, d and p, the constant Ji depends on p and r: and blows up as 
Tq — Too ^ or i?o — Roo since Ji ~ (/0~^ + . We have finished this part of the iteration since 
q = d/2 which gives in the sequel the correct exponent in the last integral. In case d = 1, 2 we have to 
observe that g = 2 so that the exponent is s = l/(po + 2(m — 1). 

This fact forces the final cylinder Qoo = (Jooj J'] x BRoo C Qq = (To,T] x B^^ to be strictly contained 
in the initial one. We obtain 



sup |f I < C] 



loc 



1 

{Rq — Ro 



+ 



1 



too — -to 



9+1 




PO+(j{m-l) 




If v^'^dxdt 




J Jqo 



PQ+ij(m-l) 



(2.37) 



Notice that Cioc only depends on m, d and p. We conclude the proof by going back from v to u, using 
the fact that u < v < u + l^hy definition of v. ^From (2.37) we easily get 

1 



sup|m| <sup|t;| < Cioc 

Qoo Qoo 

< Cioc 



9 + 1 




wP« dx dt 



pg+q{m-l) 



[[ uP«dxdt + Vol(Qc 



P0+9{™--l) 



This concludes the proof, after changing the notation, putting p = po, Roo = Ri < Ro , ?oo = Ti > To, 
and q = d/2 since we are dealing with d > 3 . □ 
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2.3 Local Smoothing effect. Proof of Theorem |2.1 



We now combine the results of Theorem 2.3 and of Theorem |2.4[ to prove the local smoothing effect 



in the form described in Theorem |2.1[ We consider u defined in (0,T) x Brq^xq), then take a smaller 
radius Ri and write Ri = {1 — e)Q and Rq = (1 + e)^: this defines p and e. We then consider the 
rescaled solution 



u{t, x) = Ku{Tt,gx + xq) 



K 



(2.38) 



with < T < r. Then, we apply to u the result of Theorem 2.4 over the cylinders Qo = (0, 1] x Bi 
and Qi = (e^, 1] x i?i-e, for some e e (0, 1) , so that the bound reads 



1^1 ^ ^loc 

sup|n| < 
Qi ^ 



£■ p+q{m-l) 



w 



Qo 



Pdxdt + Yol{Qo) 



P + 9(™-l) ^ Cloc 



9 + 1 



dx dt + LJd 



Qo 



1 

p+g(m — 1) 



since Vol((5o) = "^d • Moreover, we know that Cioc only depends on m,d,p . 

Next, we estimate the time integral, using Theorem |2.3| applied to the rescaled solution u on the balls 
Bi C -Bi+e and for times t G [0, 1] and for p > pc > 1 — m : 



/ \u{t,x)\Pdx <2^-~^ \u{0,x)\Pdx + 2^~'^[Ke,pt]^ 



(2.39) 



where 



K_ = ^voi(i?i+a^i) 



(l-m)/p ^ / 1=^=^-2 



An integration in time over (0, 1) gives 

n^o dxdt= [ [ W dxdt<2^- 



|M(0,x)|Pdx + 



2^- 



p I 1 

-Bi+e 1-m ^ 



= So[ \uiO,xWdx + ^^^ 
Jb^ e(i-m) 

and we remark that ICm,d,p and Sq only depend on m,p . 

Now we put together the above estimates, and we rescale back from u to u, also changing variable 
gx = y in the integrals, and we obtain (using K^^^ = /r) 



sup u{s,y) < K ^—3%^ 

!r2T-^lvR/. . p+q(m — 1) 



(s,?/)e{e2r,r]xS(i_^)p 



< Q 



p-\-q{m — l) ^ 

2 2±i 



cSo / |g(0,x)rdx + ^%^ + u;, 

|^z(0,2/)rdy 



p+g(m — 1) 
9l 



p+g(m — 1) 



(2.40) 



+ K2 



-"loc 



2 9+1 

^ p + g(m-l) 



m,d,p 



2p 1 
£(l-m) 



p + g(m — 1) 



In dimension d > 3 we have taken q = d/2 which allows to cancel the appearance of g in the first term 
of the right-hand side and simplify the dependence on r. We then have 



sup u{s,y)<-T^ 



B 



\uq{x)\^ dx 



Rq 



2i?„ 



+ C2 



1 

1 — m 



(2.41) 
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where we have also used (a + b)'^ < Kia'^ + K2V ■ Putting t = t we have obtained in particular the 



desired formula (2.1 ) . We last remark that e must be strictly positive even if it can be chosen arbitrarily 



small; in the limit e — > 0, the constants Ci blow up, since 

^loc , 77 ^loc 



p+q(m — 1) 



2 9 + 1 
£■ p + 9(m-l) 



and 



Co 



K2 



2 g+i 

£■ p+q(m-l) 



m,d,p 



2p 1 



p+q{m — 1) 



(2.42) 



We conclude the proof by switching to the same notations as in the statement of Theorem 2.1 just by 
substituting i?i = (1 — e)Q and i?o = (1 + £)q, it is clear that the result holds for any Ri < Rq, and 
that the constants Ci blow up when Ri Rq. We recall that Cioc only depends on m,d,p. 

We have thus proved the following result: 

Theorem 2.8 Let p > 1 if m > rric or p > pc if m < rric- Then there are positive constants Ci, C2 
such that for any < Ri < Rq we have 



sup 

(s,S/)G(io,t]x_Bflj 



B 



\uo{x)\^ dx 



+ C2 



t 



1 

1 — m 



(2.43) 



where to = [{Rq — i?i)/(2i?o)] ^ ^ and the constants Ci depend on m,d and p, Ri and Rq and blow up 
when Ri/Rq 1; an explicit formula for Ci is given by (2.42). 



The above theorem is nothing but a slightly stronger form for Theorem 2.1 we just take the limit 



proof of Theorem 2.1 □ 



i?i — > in inequality (2.43 ) to obtain (^2.lL The final expression for the constants Cj in (2.1 ) corresponds 



to the limit of Ci in (2.42) as e — > 1 (i.e. Ri — > 0) and do not depend on the radii. This concludes the 



Part III. Harnack Inequalities 



By joining together the local upper and lower estimates obtained in Parts I and II we can draw 
interesting conclusions in terms of special forms of Harnack Inequalities. These are expressions relating 
the maximum and minimum of a solution inside certain cylinders. In the standard case one has 



supu(t,x) < C infn(t, x) 



(3.1) 



see [26J and |28) . The main idea is that the formula applies for a large class of solutions and the constant 
C that enters the relation does not depend on the particular solution, but only on the data like m, d 
and the size of the cylinder R, but not on time. The cylinders in the standard case are supposed to be 
ordered in time, Qi = [ti,t2] x Br{xq), Q2 = [ti,,ti\ x Bji{xQ), with ti < t2 < ^3 < ti- 
lt is well-known that in the degenerate nonlinear elliptic or parabolic problems a plain form of the 
inequality does not hold. In the work of DiBenedetto and collaborators, see the book |17) or the recent 
work |18J, versions are obtained where some information of the solution is used to define so-called 
intrinsic sizes, like the size of the parabolic cylinder(s), that usually depends on u{tQ,XQ). They are 
called intrinsic Harnack inequalities. The authors of jl8j show that the size of a convenient cylinder 
for the Harnack inequality to hold has the form 

lR{to,XQ) = (to -cu(^o,xo)^~"i^^^o + cn(^o,xo)^"™i^^) x ^^^(xo) 
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with a fixed constant c > which depends only on m,d , that can be chosen "a priori" , but only in 
the good range nic < m < 1 . This cylinder is called intrinsic because it depends on the value of the 
solution ti at a given point {to,xo) . 

The Harnack Inequalities of |20|ll8j. in the supercritical range then read: There exist positive constants 
c and 5 depending only on m, d, such that for all {tQ^xo) S Q = (0, T) x 17 and all cylinders of the type 
IsR C Q, we have 

cu{tQ,xo)< inf u{t,x) 
xeBfi(xo) 

for all times to — Su{tQ, xq)^""™ < t < to + 6 u{to, xq)^^"^ R^. The constants 6 and c tend to zero as 
m ^ 1 or as m ^ rric ■ 

They also give a counter-example in the lower range m < rric, by producing an explicit local solution 
that does not satisfy any kind of Harnack inequality (neither of the types called intrinsic, elliptic, 
forward, backward) if one fixes "a priori" the constant c . At this point a natural question is posed: 

What form, if any, the Harnack estimate might take for m in the subcritical range < m < rUc ? 

The following is an answer this question. 

New approach. After the introduction of the lower bounds of the Aronson-Caffarelli type, it became 
clear that the size of the initial or norm in a certain ball can be used in a natural way to define 
intrinsic quantities for later times, and this is the approach the authors followed in [8] for the easier 
range rUc < m < 1 . The Harnack inequalities we derive below are based on such an idea and apply 
also for < m < rric. Indeed, if one wants to apply the result of DiBenedetto et. al. pT| UH] mentioned 
above, to a local weak solution defined on [0, T] x Q, where T is possibly the extinction time, the 
Harnack inequality of [201 HH] reads: 

There exists positive constants 5 <c depending only on m, d such that if 



cu{to,xo) < 



min{to,r - to} 
(8i?)2 



and 



dist(xo,5r2) < 



R 



(3.2) 



we then have that 



cu{to,xo) < inf u{t,x), 



for all times to — 6 u{tQ, xq)^ R'^ < t < tQ + 6 u{tQ, xq)^ R? . The constants 5 and c tend to zero as 
m ^ 1 or as m ^ rric ■ The intrinsic hypothesis (|3.2|) is guaranteed in the good range by the fact that 



solutions with initial data in L^^^ are bounded, while in the very fast diffusion range hypothesis (3.2) 
fails, and should be replaced by : 



u{t,xo) < 



d 

\u{to)\\L^(Bn)R'' 



i?2 



This local upper bound can be derived by the smoothing effect of Theorem 2.1 whenever to + ei*(io) < 



t < to + t*{to), see full details in the proof of Theorem 3.2 



The Size of Intrinsic Cylinders. We will show that the new critical time 



t*(s) 



Cm,dR^ '^^^ "'^II^(^)IIl1(Bh{xo)) 



(3.3) 
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introduced in Part I, gives the size of the intrinsic cyhnders: in the supercritical fast diffusion range 
this time can be chosen a priori just in terms of the initial datum, but in the subcritical range its 
size changes with time; roughly speaking the diffusion is so fast that the initial local information is 
not relevant after some time, which is represented by t* . We must bear in mind that a large class of 
solutions completely extinguish in finite time. We proceed next with the new results. 



Inequalities of Forward, Backward and Elliptic Type. For small times cf. Theorem 3.1 or 
for suitable intrinsic cylinders, cf. Theorem 3.2 we obtain inequalities where the infimum is taken at 
a later time than the supremum (forward Harnack inequalities), or at the same time (elliptic Harnack 
inequalities), or even at an earlier time (backward Harnack inequalities). 

Throughout this section we take < m < 1 and consider a local nonnegative weak solution u of the 
FDE defined in a cylinder Q = (0, T) x 0,, taking initial data n(0, x) = uq^x) in L^^^(J7) , with p = 1 if 
nic <m<lorp>pcifO<m< rric- We make no assumption on the boundary condition (apart from 
nonnegativity) . Also, let xq be a point in 0, and let 6-R < dist(xo, 517). As before, we let Tm be the 



so-called minimal life time, corresponding to data uq and ball Bji{xo), and we define t*(s) as in (3.3) 
and t=K = t=i,(0), which is equal or less than T^- First we prove Harnack inequalities for initial times 



Theorem 3.1 Under the above conditions, for any to € (0,t*] and < i? < min{t* — to, to/2} the 
following inequality holds 

inf u{to±e,x) >nu{to,xo) (3.4) 



where 



n =CaR- 



IkollLi(Bfl) 


m 


1 

rp {l-—m) 





I ||2P'?P 

I^OllLP(Bfl) _J_ 



t 



1 — m 





(3.5) 



and Cq depends only on m,d,p . 7i goes to zero when to 0. 



This form of Harnack estimate we propose must be called generalized, since the constant depends on 
the solution through certain norms of the data. But we remind the reader that a proper restriction of 
the class of initial data allows to control TC at any time < t < t*. 

Proof. The proof consists of two steps. 

From center to infimum. First of all we have to pass from the center to the minimum in the positivity 
estimates of Theorem 1.1. Fix a point z £ Br{xq), ad consider the following MDP, centered at z: 



dtu = A(n™) 
u(0,x) = uo{x)xbr{xo) 
u{t, x) = 



in Qt,r, = {0,T) X BgR/2{z) 

in Bn{xo), and supp(tio) C B2r{z) 

for t > and x £ dBgji/2{z) , 



(3.6) 



it is then clear that Br^xq) C B2r{z) C Bgji/2iz) C Bqr{xo). Applying the result of Theorem 
the solution u to the above minimal problem with minimal life time Tm = Tm{uo), we get 



1.6 



to 



ko 
2 



9 



-R-2R 



B2r{z) 



uq dx 



\0l{BgR/2{z)\B2R{z)) 



km,d R 



2-d(l-m) 



Uq dx 



Br(xo) 



1—m 



< T„ 



(3.7) 
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which does not depend on z{t). We then obtain: 



m — - 



uo{x)dx = ci R^'h^---Tm^ 



B2r{z) 



Br(xo) 



uq{x) dx. 



for any t G [Ojt,,,]. Since z G Br^xq) is arbitrary and does not enter neither in the above lower bounds, 
neither in the formula (3.3) for = t=i,(0), we can set z = z{t) as the point such that: 



inf u{t,x) = u(t,z{t)) > c\ 



(3.8) 



for any < t < t^:, where t^: = t*(0) is given by (3.3). Once we have obtained the result for the minimal 
Dirichlet problem we pass to a general weak solution as it has been done in Subsection 1.2 hence 
estimate (3.8) holds for any weak solution. 

Joining upper and lower estimates. Let now to G {0,t^:] and choose 9 > small so that to/2 < tQ — 6, 
and to + 9 < t^:. By the lower estimate (3.8) we know that u{tQ,xo) is positive for any to £ (0,t*] and 

(3.9) 



inf u(to ±9,x) > i?^||no|L"l,„ .T^ 

X<^Bn{xo) ^ y i H ML (Br) 

>ci2 i-'"||no|lLi(5^)r, 



m(l — m) 



(to ± 9) 



1 

1 — m 



1 2~d 1 

m(l — m) T~t 

/l ™ tn 



We now use the local smoothing effect of Theorem 2.1 for t = to G (0, t^) : 



u{to,xo) < Cs 



J. 1 — m 

''0 



i?2 



(3.10) 



We have thus proved that for to G (0,t*] and 9 < min{t* — to, to/2} we have 



I ||2pi?p 



inf u(to ±9,x) > C2 ll^o|lLi(B^)7'm™'^ '"^ R' 



x&Br{xo) 



t, 



2pi?p 
1 — m 





-1 



n(to,xo). (3.11) 



This concludes the proof, n 

By shifting the interval [0, t^:] to [to, to + t*{to)] C [0, T] we can prove a more intrinsic flavored version 
of backward-forward-elliptic Harnack inequality in the following 

Theorem 3.2 Under the above conditions, there exists constants hi,h2 depending only on m,d,p, 
such that, for any e G [0, 1] the following inequality holds 



inf M(t lb 1?, x) > /ii £1 
xeBR{xo) 



Hto)\\ 



■.Rp 



Hto)\y{BR)R'' 



(t,Xo) 



(3.12) 



for any 



to + et^{to) <t±7?<to + t,(to), t^{to) = h2R^-'^^^-"'^\\u{to)\\}- 



L^Br{xo)) 
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Proof of Theorem 3.2 Assume to = 0, the result will follow by translation in time. We continue the 
proof of Theorem 3.1 we further estimate the local smoothing effect of Theorem 2.1 for t = to G (0, t^,) : 

to 



u{to,xo) < Cs 



'i^i?I^ + l 

1 — m 



< 



\M\lp{Br) 



'to_ 



1 

1 — m 



< Ca 



et 



2pi>p 
1 — m 



1 

1 — m 



1 

1 — m 



(3.13) 



since we have put to ^ ^t* and t* = t*(0) as in (3.3). Next we use the estimate for the extinction time 
proved in Sections 1.4.1 and 1.4.2 , that can be rewritten as 

— — 2 d 

Tm''' < km,p,dR^-"' " \\uo\\lp(Br) for any p > max{pc, 1}- 



The lower estimates (3.9) becomes 



inf n(to ±0,x) > ci 



Wo\\L^{Bii)R^~ 



Ti (1 — m) 
m 



i?2 



1 

1 — m 



>C2 



\\uo\\lp(Br) 



d 

Rp 



R<^ 



i?2 



1 

1 — m 



(3.14) 



Joining now inequalities (3.13) and (3.14) we obtain that there exists constants hi , /i2 depending only 
on m, d,p, such that, for any e G [0, 1] the following inequality holds 



inf u[t ±'d,x) > /ii £1 
xeBnixo) 



u 



Rp 



|n(0)|| 



Lp(Bh) 



R'^ 



^P^P + 7^. 



{t,Xo) 



(3.15) 



for any 



et*(0) <t±i?<t*(0), u{0) = h2R 



2~d{l-m) I 



^ ^IILi(B^{xo)) 



We conclude the proof by translating the result from [0,t=K(0)] to [toi^o + ^*(*o)]) that we know to be 
included in [0, T] as explained in Part I . □ 



Remarks, (i) Estimate (3.12) is completely of local type, since it involves only local quantities. In 
the supercritical range rric < m < 1, we can let p = 1 in (3.12) to get 

inf n(t lb -i?, x) > hie ^-"^u[t, xq) 
xGBn{xo) 

for any to + £t^:{to) < t ± < to + t*(to)- In this way we recover the above mentioned results of 
DiBenedetto et al., cf. |20| I18j. This theorem complements and supports the lower Harnack inequality 
(1.9) of Part I, in which the constant does not depend on uq. 

Joining the upper and lower estimates for the Cauchy Problem, we obtain the Global Harnack principle 
as the authors did in [8] . 

(ii) In the subcritical range < m < rric, the Harnack estimates cannot have a universal constant 
independent of uq, as already mentioned, cf. also in [18] for a counterexample. We have proved that 
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if one allows the constant to depend on the initial data, then it is possible to obtain intrinsic Harnack 



inequalities, and the price we pay is having the minimal life time in the constant, as in Theorem 3.1 



but this information is a bit unpractical and we replace it with some local L^-norm with p > pc of the 
initial datum. 

(iii) We have shown that the size of the intrinsic cylinders is always proportional to a ratio of local 
norms. Note that in the supercritical range it simplifies and only depends on the local norm. 

(iv) The quantity e represents an arbitrary small waiting time, that is needed in order for the regular- 
ization to take place and to allow quantitative intrinsic Harnack inequalities. 

(v) Backward Harnack inequalities are a bit surprising, but they reflect a typical feature of the fast 
diffusion processes, that is the extinction phenomena, namely 



inf u[t — -djx) > /ii £1 

xi^Bji{xo) 



k(io)|| 



Rv 



Hto)\\ 



Lf(i?ij) 



[t,Xo) 



(3.16) 



for any 

to + eU{to) <t-^ <to + U{to), t*(to) = /l2i^^"'^^^""'^||n(^o)|lLT(B^(^o))• 
This inequality is compatible with the fact that the solution extinguish at some later time, remaining 
strictly positive before. This backward inequality is typical of singular equation and can not hold for 
the degenerate -porous media- case m > 1, neither for the linear heat equation case, m = 1. 

The same remark applies for the Elliptic Harnack inequality, that is when -i? = 0. 

An alternative form of Harnack Inequalities. We provide a form of Harnack inequalities of 
forward, backward and elliptic type, avoiding the intrinsic framework, and the waiting time £ £ [0, 1]. 

Theorem 3.3 Under the above conditions, there exists positive constants Ci, C2 and /12 depending only 
on m, d and p such that 



Ci 



sup u{t,x) < \\u{tQ)\\l^^'' +C2 



Hto)\\ 



Hto)\\ 



d 

Rp 



inf u(t±^,x) (3.17) 

x£Br 



for any 



0<to<t±^ <to + t^{to) <T, t^{to) = h2 R^-'^^^-"'^\uito)\\]j 
where -dp = l/(2p — d{l — m)). 



m 

LHBh(xo)) 



Proof. First we observe that we can pass from the center xq to the supremum in the upper estimate of 



Theorem 2.1 by doubling the radius of the ball on the right hand side, namely 

There exist positive constants Ci, €'2 depending only on m, d, such that for any t,R > we have 



sup lift, x) < -^^llufto)!!?^'!?? / ^^ + Co 



i?2 
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Joining the above inequality with the lower bound of Theorem in the form of (3.14), we obtain the 
inequality (3.17) for to = 0. We conclude the proof by shifting the interval [0,t*(0)] to [toi^o + i*(io)] j 
that we know to be included in [0, T] as explained in Part I . n 



Remark. In the good fast diffusion range we can let p = 1, so that inequality (3.17) reads 

Ci 



sup u{t, x) < 
x&Br 



td^, Ht0)\\f+C2jg^u{t±i},x). 



4 Concluding section 

First, we sketch a panorama of the local estimates in the different exponent ranges m < 1 and for 
different integrability exponents p > 1 that may serve as further orientation for the reader. Then we 
make a series of general comments, and finally we review related works. 



4.1 Panorama and open problems 

The values of rric and pc are defined in the Introduction. 




(I) 



Good Fast Diffusion Range: m G {rUc, 1) and p > 1. Here, the local smoothing effect holds, cf. 
Theorem |2.1[ as well as the Reverse smoothing effect and the global Harnack principle, proved in 
[8] , see also [9] . As for the result of the present paper, we have provided a different proof of the 



positivity result in Theorem 1.3 We have also proved intrinsic Harnack inequalities of forward. 



(11) 



elliptic, or backward type, cf. Theorems 3.1 3.2 and 3.3 recovering the existing results. Finally, 
for times close to the extinction time, in case extinction occurs, the authors show in |10| |9]. that 
Elliptic Harnack inequalities hold up to extinction. 

Very Fast Diffusion Range: m £ (0, nic) and p > Pc > 1 • The local smoothing effect of Theorem 
|2.1| holds, as well as the lower estimates of Theorem |1.1[ These are the only known local positivity 
and smoothing results in this range. For any positive time we have the Aronson-Caffarelli type 
estimates of Theorem |l.l[ We have also proved intrinsic Harnack inequalities of forward, elliptic. 



or backward type, cf. Theorems 3.1 , 3.2 and |3.3[ These are the only parabolic Harnack inequalities 
known in this whole range. 
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An open problem is to pass from local to global estimates in this very fast diffusion range. For 
m > mc this is done in [8] in the form of a global Harnack principle. 

Another open problem is to find the rate of convergence to an appropriate extinction proGle for 
the Dirichlet problem on domains. A minimal rate of extinction can be obtained by our intrinsic 



Harnack estimates, of Theorem 3.2 details will appear separately. 



(c) Critical case: m = mc and p > pc = 1- The local upper and lower estimates of zone (II) apply, 
as well as the Harnack inequalities. In an upcoming paper we will show how to pass from local 
to global estimates, obtaining global lower bounds with super-exponential time decay. 

(Ill) Negative exponent range: m <0 with p > pc- No positivity result is known in this range, and the 
technique developed in this paper does not allow to treat this case. Recall that solutions of the 
homogeneous Dirichlet problem on bounded domains vanish instantaneously (hence, there is no 



actual solution). The local upper bound of Theorem 2.1 is still valid, and this is the only known 
local upper estimate in this range. It is an open problem in this case to End positivity and a 
posteriori Harnack inequalities, if any. 

(IV) Negative range: m <0 with p < pc- The smoothing effect is not true, since initial data are not in 
with p > Pc, cf. [31], and the solutions of the Cauchy problem with data in L^(]R'^) will vanish 
instantaneously, setting a strong limitation to positivity results, which must be based of course 
on local bounds. Again, solutions to the homogeneous Dirichlet problem on bounded domains 
vanish instantaneously. It is an open problem to find positivity estimates in this range, if any. 
In general, Harnack inequalities are not possible in this range since solution may not be (neither 
locally) bounded. 

(V) Very Fast Diffusion Range, m £ (0,mc) with small integrability exponent p G [l,Pc]- It is well 
known that the smoothing effect is not true in general, since initial data are not in with 
p > Pc, cf. [H [T2l |3T] . The lower estimates in this case are given by Theorem These are the 
only known local positivity results in this whole range. In general, Harnack inequalities are not 
possible in this case since solution may not be (neither locally) bounded. 

4.2 Some general remarks 

• We stress the fact that our results are completely local, and they apply to any kind of initial- 
boundary value problem, in any Euclidean domain: Dirichlet, Neumann, Cauchy, or problem for 
large solutions, namely when u = +oo on the boundary, etc. Natural extensions are fast diffusion 
problems with variable coefficients and fast diffusion problems on manifolds. 

• We calculate (almost) explicitly all the constants, through all the paper. 

• Our positivity and Harnack inequalities generalize the results of \19\ [20l [T8] , valid only in the 
good fast diffusion range, in the sense that we recover their result with a different proof, and we 
extend it quantitatively to the very fast diffusion range. 

• More specific information can be obtained when we restrict our attention to particular classes 
of solutions, like the solutions of the Cauchy Problem, the solutions of the initial and boundary 
value problem on a bounded domain with n = on the boundary, or with u = -|-oo on the 
boundary. All these solutions have additional properties that can be exploited. We will not enter 
into those topics for reasons of space. 
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• We will not enter either into the derivation of Holder continuity and further regularity from the 
Harnack inequalities. This is a subject extensively treated in the works of DiBenedetto et al., see 
|2H \T7\ [T8] and references therein. 

• Actually, the fast diffusion equation can be well-posed even for Borel measures (i.e., not locally 
finite measures) as initial data. This is done for the Cauchy problem in the whole space in 
|12] . but in that case the smoothness of the constructed solutions is lost and the concept of 
extended continuous solution has to be introduced. The problem turns out to be well-posed in 
that extended class in the good range nic < m < 1. 

• The main ideas of this paper can be extended to related equations, such as the p-Laplacian, but 
the technical details may not be immediate. 



4.3 Elliptic connection 

Part of the difficulties of the FDE in the lower range of m can be explained by the intimate relation 
of the equation with semilinear elliptic theory. Indeed, if we try the separation of variables ansatz 
u{t,x) = H{t)F{x) on equation (0.1) we easily get the possible formulas Hi{t) = t^/(^~"^) or H2{t) = 



(T — t)^^^^ and then G = F"^ satisfies the elliptic equation 

-AG±cG^ = 0, withg = — , (4.1) 

m 

so that g>lifO<m<l. The constant is c = m/(l — m), and the sign it corresponds to the choices 



Hi or H2 respectively. It is well-known that the theory of equation (4.1) is difficult for large values 
of q, notably for q > qs = {d + 2)/{d — 2). The exponent corresponding to uIq is qc — d/(^d — 2), a 
lower exponent that appears sometimes in the study of singularities. Note that when m is negative G 
is an inverse power of F; moreover, q is negative]^ The FDE-elliptic connection extends to the study of 
self-similar solutions of different types, which is a fundamental tool of the FDE theory and is described 
in [25] and [31 J among other references. 

It is interesting to check the correspondence of our time evolving results with the theory of elliptic 
equations. An easy way of doing that it to apply the results to special solutions. The simplest case, 
i.e., stationary solutions, is too simple, hence we prefer to try the separate variable solutions 

u{t,x) = (T-t)i/(i-™)F(a;). 



Putting U = F"' and q = 1/m we get the elliptic equation AU + cW^ = 0, as in (4.1 ). It can also be 
written as 

AU + a{x)U = 0, a{x) = cU^-^ = cF^-"^ . 



Let us check the upper estimate ( |2.1[ ) of Theorem 2.1 Using the separate variable form of u we 
immediately see that the time dependencies disappear (a confirmation of the correct scaling of the 
formula) and the we obtain a local boundedness result for U of the form 

U{x) < C||?7|p^) + C72i?-2P/(1-P)^ 

on the condition that F G L^^^(r2) with p > Pc which means that a(x) G L[^^(J7) with r > d/2, a 
classical condition. This is for us another way of checking that pc is a natural exponent. 



■^We leave to the reader the exponential formulas that are obtained for m - 
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A similar conclusion can be derived from application of the lower estimate (1.1 ) of Theorem |1.1[ We 
leave the details to the reader. The elliptic conclusions can also be checked on selfsimilar solutions of 
the type u{t, x) = {T - t)°F{x (T - tf), as the ones considered in [3l] . 



4.4 Short review of related works 

The range m < rric has remained outside of most of the publications on the questions of positivity 
and Harnack estimates. For positivity and boundedness, let us first mention the works of Bertsch and 
collaborators [5j 6^ who treat the equation satisfied by the so-called pressure variable v = c/u^~^, i.e., 
an inverse power of u. It covers the equivalent to the whole fast diffusion range m < 1 in terms of 
viscosity solutions; the questions are somewhat different from our program. 

We list next and comment on a number of previous results on the subject of Harnack inequalities for 
the fast diffusion equation. 

• DiBenedetto and Kwong proved in fl9] that in the good fast diffusion range rUc < m < 1 intrinsic 
Harnack Inequalities of forward type do hold, under the positivity assumption that there exists 
a point (to,xo) such that u{to,xo) > 0. This value controls from below the infimum in a small 
ball at a later time, with sizes depending on u(xo, to). 

• Later on DiBenedetto, Kwong and Vespri [2^ improved on the previous result proving the Global 
Harnack Principle in the wider range = < m < 1, by means of comparison with the 
separation of variable solution, always under the assumption of positive solutions and a stronger 
assumption on the initial data, namely G Wq'^(O). These estimates are global in space but 
not in time since the constants blow up as t — > 0. Hence the interest in combining them with 
information we provide for all small times in direct dependence of the local norms of the initial 
data. 

• More recently, DiBenedetto, Gianazza and Vespri, |18], extended the results to the variable 
coefficient case in the form Harnack inequalities which are of Forward, Elliptic and Backward 
type; it applies in the good FDE range, and always under the positivity assumption for some 
(xo, to)- Some of these estimates had been proved by the authors in the constant coefficient case 
in [8] , see also [9]. 

• The above mentioned Harnack inequalities imply Holder continuity of the solution and sometimes 
analyticity, cf. [201 [T8]. 

• The power nis = {d — 2)/{d + 2), has been studied by Del Pino and Saez [15], as part of the study 
of the asymptotics of the evolutionary Yamabe problem. They perform the transformation into 
a fast diffusion problem posed on the sphere via stereographic projection, which is possible for 
this exponent. They get an elliptic Harnack inequality which holds for a good class of solutions, 
but they do not prove a parabolic Harnack inequality. 

• None of the above quoted papers considers the problem of finding Harnack inequalities when the 
time approaches the finite extinction time (if there is one). This has been done by the authors 
in [10| |9] , showing that Elliptic Harnack inequalities hold up to the extinction time. The proof is 
completely different, we draw fine asymptotic properties, by a careful analysis of the extinction 
profile. 
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Summing up, two results are known in the lower range m < me- |15] that applies for m = rng, and 
[20j that applies for < m < 1. They both refer to a different point of view. 



Appendixes 



Al Aleksandrov's Reflection Principle 

Here we state the Reflection Principle of Aleksandrov in a slightly different form, more useful to our 
purposes. We already used this proposition, in [8]. Other forms of the same principle, in different 
settings can be found, for example in [22]), Proposition 2.24 (pg. 51) or in P], Lemma 2.2. 



Proposition 5.1 (Local Aleksandrov's Reflection Principle, [8j ) Let Bxr„{xo) C M*^ be an open 
ball with center in xq G M'^ of radius XRq with Rq > and A > 2. Let u be a solution to problem 

' dtu = A(n"^) in (0, +cx)) x BxR,ixo) 

< u{0,x) = uo{x) in BxRoixo) (5.2) 
u{t,x) =0 for t > and x E dBxR^{xQ) 

with supp(tio) C Br^{xo). Then, for any t > one has: 

u{t,xo) > u{t,X2) 

for anyt>0 and for any X2 G Ax,r^{xo) = BxRoixo) \ B2Ro{xo). Hence, 

> |A,^(.o)r' / .(t..)d.= / „(*,.)d. (5.3) 



The proof of this result can be found in the Appendix of ^ . 



A2 On the extinction time 



The extinction time plays a role in our study of positivity in Part I. We devote some paragraphs to 
comment on its occurrence in FDE. It has been observed in the literature, cf. [3 [U [31], that Lemma 



2.2 can be used to obtain lower bounds on the finite extinction time T = T{uq). Indeed, just let t = 0, 



s = T and w = 0, in (2.2), to get 

ci(m,d,A)i?'^(i— )-2||no||i7('^^)< T (5.4) 

Suppose that = M*^. When rric < m < 1, it is easy to see that 2 — d{\ — m) > 0, so that, simply by 
letting R — > oo, we see that T{uo) — > +oo if uq G L^(M'^), that means that solutions corresponding to 
initial data uq G L^(M'^) do not extinguish in finite time, i.e., there is global positivity. 

Such a situation does not occur for solutions defined in all of M'^ when < m < rric. Indeed, in that 
fast diffusion range solutions with initial data in L-'^(M'^) may extinguish. The question is studied in 
Chapters 5-7 of |31| where upper and lower estimates for T in terms of the data are obtained for the 
problem posed in the whole space with m < rric- 
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On the other hand, solutions extinguish in finite time for the Cauchy-Dirichlet problem posed in a 



bounded domain with zero boundary conditions for all < m < 1, cf. subsections 1.4.1 1.4.2 or [?. I16|. 



An estimate similar to (5.4) cannot be valid for m < since it is known that the Cauchy problem 
does not admit solutions with data no S L^(M'^) |30], or even in LP(M"') with p < pc [13] for m < 0. We 
may say in these cases that T = 0. Actually, when one tries to solve approximate problems with, say, 
strictly positive and bounded boundary data and pass to the limit, the approximate solutions converge 
to zero uniformly in cylinders of the form Q-^ = (t, oo) x M'^. This can be summed up as the formation 
of an initial discontinuity layer, cf. |3T]. No solutions exist either for the Cauchy-Dirichlet problem 
posed in a bounded domain with zero boundary conditions when m <0. There is also a peculiar case, 
m = and d = 2, where there exist solutions but the waiting time can be fixed a priori independently 
of the initial data [31], hence no estimate as above is possible either. 



Note finally that the estimate of Theorem 2.3 for p > 1 cannot be used to obtain lower estimates 
on FET, since they hold only for t > s >0. 



A3 Details of the iterative calculations of Subsection 2.2, Step 3 



We show here in detail of a calculation used to pass to the limit when A; — > cxd in the inequality (2.34) 



in the proof of Theorem |2.4[ We adopt the notations used there. 

k 



j Pk+i < exp 



(1 + l)k-3 



Pk+i 



< exp 



mj 



1 + 



fc+i 



+ (g + 1)(1 - m) 



< So exp 



< si exp 



k 

2^(l+^)-^~Mog(i) 

oo ^ 



< si exp 



k 



2(9+1) 



where Sj, i = 0, 1 are positive numerical constant. □ 
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